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UNDERWATER  SOUND  TRANSMISSION 


Preface 


The  subject  "Underwater  Sound"  covers  phenomena  as  complex 
and  diverse  as  the  imagination  and  patience  of  the  investigator  will 
allow.  An  important  part  of  the  subject  deals  with  the  sonar  parame¬ 
ters,  in  the  form  of  loss  or  strength  coefficients.  During  the  twenty 
years  that  appreciable  effort  has  been  expended  in  experimentation 
and  analysis,  the  great  emphasis  has  been  in  two  areas. 

a)  Quantitative  investigation  of  the  sonar  parameters. 

b)  Quantitative  investigation  of  the  ocean  environment. 

Knowledge  of  the  environment  has  been  applied  to  the  estimation  of 
the  parameters,  and  to  the  description  of  other  aspects  of  underwater 
sound. 

The  principle  effects  of  the  environment  which  determine  losses 
in  a  propagating  sound  wave  were  identified  by  the  early  part  of  World 
War  IL  The  work  since  that  time  has  served  to  provide  an  increased 
understanding  of  the  physical  processes  involved,  and  to  make  possible 
more  detailed  and  reliable  calculations  of  transmission  loss. 

The  first  concern  of  the  user  or  the  designer  of  underwater  sound 
equipment  is  with  energy,  or  power.  These  things  establish  the  general 
operational  features,  size  and  weight  of  equipment,  and  class  of  compo¬ 
nents.  The  state  of  present  knowledge  permits  the  evaluation  with  reason¬ 
able  accuracy,  of  energy/power  levels  under  prescribed  conditions.  Un¬ 
fortunately,  conditions  prescribed  must  be  quite  exact,  if  definite  energy 
calculations  are  required.  This  is  true  because  present  descriptions  of 
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the  sound  field  are  synthetic,  in  the  sense  that  the  whole  picture 
is  made  up  of  many  parts,  each  part  contributing  an  important 
element  to  the  calculation.  Each  part  also  contributes  its  share 
to  the  uncertainties  of  the  reliability  of  the  calculation. 

Table  I  outlines  the  conditions  which  must  be  specified  for 
a  calculation  of  transmission  loss,  as  an  example.  It  is  clear 
that  a  specific  calculation  will  apply  only  to  the  specific  conditions 
chosen.  Now,  of  course,  no  one  has  more  than  passing  interest  in 
such  a  specific  situation.  Therefore,  in  practice  one  either  makes 
many  calculations,  corresponding  to  many  conditions,  and  then 
averages  them  in  some  sense,  or  else  chooses  a  "representative" 
condition,  and_  regards  it  as  typical. 

The  process  of  research  and  application  as  discussed  above 
amounts  to  measurement  of  the  sound  field,  its  analysis  into  many 
constituents,  and  then  a  synthesis  of  the  predicted  field  by  a  sort  of 
reverse  process.  But  the  reverse  does  not  lead  back  to  the  original 
measurement,  even  in  principle.  All  of  the  detail,  and  most  of  the 
features  of  the  sound  field  have  been  irreversibly  discarded.  There¬ 
fore,  no  answers  can  be  provided  to  questions  concerning  other  fea¬ 
tures  of  the  sound  field,  and  if  the  questions  must  be  answered,  the 
process  of  research  and  application  must  be  renewed  and  redirected. 

Today,  the  questioners  are  pounding  on  the  door,  seeking  infor¬ 
mation  vital  to  the  problems  of  present  operational  performance  and 
future  design.  They  are  concerned  for  the  most  part  with  features  of 
the  sound  field  which  are  simply  not  contained  in  the  sonar  parameters. 
There  are  two  ways  in  which  answers  to  such  questions  can  be  sought 
at  this  time.  One  is  the  direct  way,  going  to  sea  again,  with  new  kinds 
of  instrumentation,  making  new  analyses  and  building  up  a  new  picture 
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of  the  kind  that  seems  important.  The  second  way  is  synthetic, 
using  a  physical  description  of  the  ocean  as  the  basis  for  the  neces¬ 
sary  equations  to  describe  the  required  features  of  the  sound  field. 
Both  ways  are  useful;  both  are  being  employed  today. 


Unfortunately,  the  kind  of  experiments  being  conducted  today 
cannot  be  regarded  as  satisfactory.  On  the  "synthetic"  side,  con¬ 
sidering  the  very  large  expenditures  of  effort  in  the  past,  future  re¬ 
sults  cannot  be  expected  to  provide  more  than  marginal  increases  in 
knowledge,  except  in  special  cases.  On  the  "direct"  side,  there  is 
the  inescapable  dilemma  that  vastly  increased  detail  is  required,  on 
a  subject  which  is  already  almost  unmanageable  because  of  detail. 
The  fact  is,  new  ideas  and  concepts  for  characterizing  the  sound 
field  are  sorely  needed.  What  are  the  essential  features  of  the 
sound  field,  ich  both  must  and  can  be  determined  ? 

This  > o port  pp»e out physical  representation  of  the  ocean, 
together  with  constants  and  equations  sufficient  to  permit  a  reason¬ 
able  calculation  of  the  average  sound  field.  In  addition,  it  discusses 
other  features  of  the  sound  field,  and  hints  at  some  new  methods  of 
calculation  using  the  physical  picture. 


>  b«H«f  reviev^of  the  state  of  knowledge  con¬ 
cerning  underwater  sound  transmission  from  its  status  in  1946,  to  its 
present  state.  A  discussion  of  recent  advances  and  studies  in  several 
important  areas  is  made.  An  oceanographic  acoustic  description  of 
the  physical  environment  of  underwater  sound  is  proposed  in  outline 
form. 
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In  Chapter  2,  the  proposed  description  is  developed  quantita¬ 
tively,  and  detailed  methods  are  given  for  computing  the  average 
sound  field* 

Chapter  3  discusses  aspects  of  the  sound  field  which  are  not 
covered  by  the  quantitative  model,  either  in  principle  or  because 
of  practical  difficulties.  These  include  mean  variability  of  the 
sound  field,  fluctuations  in  the  sound  field,  turbulence  effects  and 
other  measures  of  the  30\md  field. 

In  Chapter  4,  recommendations  for  the  direction  of  future  re¬ 
search  in  underwater  sound  propagation  are  given. 

To  the  reader  not  concerned  with  detail,  it  is  suggested  that 
Chapters  1  and  4  be  perused.  . 
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Chapter  1 


INTRODUCTION 


1.  1  Orientation 

Since  the  end  of  the  war,  only  two  status  reports  on  the  sub¬ 
ject  of  underwater  sound  transmission  have  been  published.  These 
include  the  NDRC*  reports  of  1946,  and  the  Geological  Memoir  No. 

27^  by  Ev/ing,  Worzel,  and  Pekeris  of  1948.  The  present  report 
reviews  the  progress  made  in  the  field  since  that  time. 

Rapid  strides  have  been  made  in  the  understanding  and  pre¬ 
diction  of  transmission  loss.  The  previous  sunamaries  consisted 
of  collections  of  articles  by  various  authors  with  diverse  points  of 
view,  Unfortvinately,  gaps  existed  in  available  data  and  in  the  gener¬ 
al  imity  of  the  subject  matter,  Mfhile  much  remains  to  be  done  in 
both  experiment  and  concept,  we  have  strived  to  achieve  tmity  and 
completeness  in  treatment  and  to  establish  an  acoustical  description 
of  the  ocean  into  which  most  of  the  present  and  future  data  can  be 
placed  without  undue  strain  on  the  structure. 

The  present  state  of  agreement  between  the  "oceanographer's 

ocean"  and  the  "acoustician's  ocean"  has  allowed  the  formulation  of 

an  acoustic  description  of  the  ocean  as  a  framework  for  additional 

o 

experimental  design  and  operational  applications.  The  cells  of  the 
descriptive  model  must  be  fine-grained  enough  to  admit  data  from 
general  research  and  yet  large  enough  to  permit  easy  assembly  and 
integration  for  specific  operational  applications  which  exist  now  and 
those  which  will  be  required  in  the  future.  The  proposed  description 
of  this  model  is  shown  in  Table  I, 


Table  I 


1.  0 
1.  1 
1.  2 
1.  2.  1 
1.  3 
1.  4 

1.  5 

2.  0 

2.  1 
2.  1.  1 
2.  1.  2 
2.  1.  3 
2.  1.  4 
2.  2 
2.  3 

2.  3.  1 

3.  0 

3.  1 
3.  2 

3.  3 

4.  0 
4.  1 

4.  1.  1 
4.  1,  2 
4.  2 
4,  2.  1 

4.  2.  2 

5.  0 

5.  1 
5.  1.  1 
5.  1.  2 
5.  2 

Koto: 


PHYSICAL  ENVIRONMENT  OF  UNTJER WATER  SOUND 

e 

Outline 


Regions 

Atmosphere 
Ocean  Surface 
Ice  cover 
Hydrosphere 
Ocean  floor 
Ocean  sub-bottom 

Geophysical  Properties  -  (point  functions  of  space-time 

in  each  region) 

Physical  properties 
Depth 

Velocity  (wind  or  current) 

Temperature 

Pressure 

Chemical  properties  ( composition) 

Biological  properties 

Presence  of  species 

Acoustic  Properties  -  (point  functions  of  the  geophysical 

properties  and  the  acoustic  parameters) 

Phase  velocity 
Absorption  coeffigient  ^ 

Boundary  reflection  and  scattering  coefficients 

Parameters  (independent  variables) 

Acoustic 

F  requency 
Wave  amplitude 

Sona.r 

Source  location 
Receiver  location 

Sound  Properties  -  (point  functions  of  the  acoustic  properties 

and  sonar  parameters) 

Sound  pressure 

Mean  value 
Mean  Variability- 
Transmission  Time 

The  model  may  be  extended  to  other  sound  properties  such  as  rever¬ 
beration.  etc.  ,  by  adding  to  items  4.0  (possibly)  and  5.0,  above. 


1.  2  Review  of  Present  Status 

1.  2.  1  General  Comments 

The  chief  advances  since  1946  consist  of  the  probing  of  the 
sound  field  as  a  function  of  near- surface  depth  for  both  source  and 
receiver;  a  better  recognition  of  the  role  of  surface  scattering  in 
the  attenuation  of  sound  in  surface  sound  channels;  establishment 
of  the  forrr^  of  the  frequency  and  temperature  dependence  of  absorp¬ 
tion  and  its  relation  to  the  MgS04  ions;  the  establishment  of  prelimi¬ 
nary  curves  for  estimating  bottom  loss  as  a  function  of  grazing  angle 
and  frequency;  some  estimate  of  the  gains  and  geometry  of  conver¬ 
gence  zone  fields;  a  summary  of  shallow  water  propagation  results 
in  terms  cf  acoustic  frequency,  water  depth,  and  bottom  type;  studies 
concerning  the  relationship  between  thermal  microstructure  and  sound 
fluctuations;  and  the  direct  measurement  in  the  laboratory  of  the  speed 
of  sound  in  sea  water,  correcting  Kuwahara' s  surface  values, 

A  great  deal  of  measurement  work  continues  to  be  done,  but  the 

emphasis  is  still  on  the  specific  needs  of  operational  equipment.  In- 

o 

formation  about  variability  in  both  time  and  space  has  become  even 
more  important  and  the  three-dimensional  representation  of  the  oceans 

is  needed  to  describe  it  as  a  medium  for  practical  acoustic  transmission. 

o 

1,  2.  2  Speed  of  Sound  and  Absorption  Measurements 

With  the  development  of  the  Greenspan-T schiegg^  sound  velocity 
meter  into  a  field  instrument,  increasing  effort  is  being  given  to  the 
actual  measurement  of  the  speed  of  sound  and  its  fluctuations.  It  has 


been  used  in  connection  with  propagation  measurements  at  seau  While 
controversy  still  exists  regarding  the  exact  nature  of  the  sound  speed 
which  this  device  measures,  the  formulas  of  W.  Wilson'*"”  have  been 
accepted  for  the  conversion  of  oceanographic  measurements  of  temper¬ 
ature  and  salinity  as  a  function  of  pressure.  These  expressions  were 

determined  from  laboratory  measurements,  but  have  been  fairly  well 

« 

verified  in  the  field.  Accurate  velocity  measurements  are  important 
for  determining  range  to  convergence  zones,  in  long-range  communi¬ 
cations  and  in  fathometry.  Mackenzie^  has  furnished  a  gravity  correc¬ 
tion  for  the  density  of  sea  water  in  converting  pressure  to  ocean  depth, 
Maunsell  and  Scrimger®  have  shown  that  Kuwahara's  formula  yields 
ranges  too  short  for  the  observed  convergence  zones  in  the  North-East 
Pacific, 

Since  absorption  is  an  important  part  of  attenuation  in  sea  water, 
a  more  precise  determination  of  that  quantity  will  aid  in  ui^er standing 
^he  residual  term,-  or  the  attenuation  due  to  other  factors.  Much  work 
has  been  done  relating  the  absorption  properties  of  sea  water  to  those 
of  magnesium  sulphate  in  solutiorj.  While  there  is  reasonable  agree¬ 
ment  with  regard  to  the  data  on  relaxation  frequencies  and  the  form  of 
the  expression  and  consequent  temperatflre  dependence,  the  value  for 

the  temperature  independent  coefficient  appears  to  be  in  dispute.  There 

o 

is  a  35%  difference  in  determinations  of  this  coefficient  which  should  be 
resolved.  There  are  only  two  large  sets  of  measurements  at  sea  from 
which  absorption  characteristics  have  been  determined.  The  first  set, 
that  of  Marsh  and  Schulkin,’'  is  in  moderate  agreement  with  that  of 
Del  Grosso, “  and  of  Kurtze  and  Tamn^^.’*^  Measurements  for  the  second 
set,  that  of  Murphy  and  Garrison,*^  are  at  higher  frequencies,  and  were 
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made  in  the  low  salinity  water  of  Puget  Sound.  Some  suspicion 
exists  here,  because  even  though  they  agree  with  Beyer's  nomo¬ 
gram.  Beyer^acknowledges  that  the  consistency  of  the  different 
salts  in  sea  water  fails  in  water  of  low  salinity,  and  absorption 
has  been  found  to  be  proportional  to  salinity.  Thus  a  study  of  the 
temperature  independent  term  in  sea  water  as  opposed  to  mag¬ 
nesium  sulphate  solution  can  turn  out  to  be  quite  important. 

Del  Grosso^^  has  computed  a  maximum  dispersion  of  0.024 
m/sec  from  an  expression  by  Kneser  that  the  percentage  disper¬ 
sion  is  proportional  to  the  maximum  absorption  at  the  relaxation 
frequency.  The  relationship  between  the  temperature  pressure, 
and  salinity  dependence  of  sound  speed  and  absorption  has  not  been 
worked  out  yet.  Nor  has  the  pressure  dependence  of  absorption  in 
sea  water  been  measured,  although  this  has  been  done  for  pure  water, 

1.  2,  2,  1  Attenuation  Coefficient 


One  of  the  earliest  techniques  of  describing  propagation  loss 
over  a  path  was  to  fit  data  to  an  equation  which  had  a  term  considering 
spreadii^  loss  and  a  term  in  which  the  loss  was  proportional  to  range. 

The  coefficient  of  this  latter  term  is  called  the  attenuation  coefficient 

e 

and  is  frequency  dependent.  Two  expressions  for  the  attenuation  co¬ 
efficient  which  are  in  fairly  widespread  use  are  those  of  Horton'*’ and 
Sheehy  and  Halley*^  .  The  Horton  formula  fits  empirical  data  above 
5  kc/s  to  the  same  frequency  dependence  as  is  given  by  the  absorption 
formula,  Sheehy  and  Halley  found  a  ^/z  power  law  dependence  on  fre- 
quency  from  40  cps  to  lOO  kc/s.  Recently  E,  L,  Peterson  developed 
an  equivalent  network  for  sea  water  using  the  Horton  formula.  He  showed 
that  the  impulse  response  for  the  network  agreed  with  the  properties  of 
explosion  impulses  in  the  seau 


1.  2.  3 


Ocean  Boundaries 


1.  2.  3.  1  Surface  Scattering 

Surface  scattering  leads  to  a  large  attenuation  term  when 
coupled  with  the  predominant  occurrence  of  surface  sound  channels. 
Marsh  et  al^'^  have  presented  a  solution  to  these  problems,  in¬ 
cluding  application  of  the  Neumann-Pier son  expression  for  ocean 
wave  spectra.  The  formulation  also  allows  the  determination  of 
the  sound  intensity  fluctuation  distribution.  Marsh^^’^^  has  also 
applied  his  analysis  to  surface  reverberation  and  has  shown  how  to 
deduce  ocean  wave  spectra  from  shot  spectra.  There  is  evidence 
pointing  to  the  existence  of  other  factors  which  must  be  considered 
to  obtain  a  better  understanding  of  the  nature  of  the  propagation 
mechanism.  For  example,  Urick^'*  has  observed  scattering  returns 
emanating  from  a  point  well  below  the  depth  of  the  waves  on  the  sea 
surface.  In  addition,  Urick  and  Hoover^®  have  observed  an  indepen¬ 
dence  of  back- scattering  intensity  below  a  10®  grazing  angle  in  rela¬ 
tively  rough  water  which  they  attribute  possibly  to  bubbles.  Garri¬ 
son,  Murphy  and  Potter^  observed  an  extremely  sharp  falloff  of 
back- scattering  at  these  low  angles  when  the  water  surface  was  com¬ 
paratively  undisturbed, 

1.  2.  3.  Z  Under  Ice  Propagation 

Several  reports  on  transmission  experiments  under  ice  have 
appeared  recently,  and  the  subject  is  receiving  accelerated  attention. 
Work  of  the  group  at  Lamont  Geological  Observatory^^  and  Naval  Re¬ 
search  ifstablishment  ( Canada)^®  has  served  to  demonstrate  several 
important  features  of  the  under  ice  situation,  foremost  among  them 


being  the  sharp  cutoff  in  propagation  above  50  or  100  cps.  This  has 
been  attributed  to  the  scattering  effects  of  the  rough  underside  of  the 
ice.  Coupled  with  this  effect  is  the  loss  of  high  order  modes,  when 
the  refraction  is  upward,  as  it  usually  is.  The  resxilt  is  a  situation 
which  can  be  characterized  by  a  few  modes  giving  rise  to  the  effects 
of  dispersion  and  modal  interference. 

On  the  quantitative  side,  the  scattering  losses  under  ice  can  be 
computed  from  reasonable  assumptions  about  the  under-ice  roughness, 
using  formulas  developed  for  sea  surface  scattering.  Mellen^’  has 
made  extensive  calculations  interpreting  the  mode  patterns. 

1.  2.  3,  3  Acoustic  Properties  of  Ocean  Bottom 

The  effects  of  the  bottom  in  its  influence  on  the  sound  field  vari¬ 
ability  are  being  studied  by  various  groups.  An  understanding  of  the 
acoustic  properties  of  the  ocean  bottom  has  become  quite  important. 
Attempts  are  being  made  to  separate  the  various  factors  affecting 
apparent  bottom  losses  such  as  sound  velocity  structure,  pulse  widths, 
and  scattering  and  absorption,  in  addition  to  the  varying  nature  of  the 
bottom  itself.  Fry  and  Raitt*®  have  started  a  classification  of  bottom 
types  into  "fast"  and  "slow".  A  dependable  estimate  of  avera^  bottom 
loss  as  a  function  of  frequency  and  grazing  angle  is  required,  as  well  as 
a  measure  of  the  dispersive  and  delay  properties. 

Hamilton^’”^'*  et  al  at  NEL  have  made  laboratory  measurements  of 

bottom  properties  by  bringing  back  samples  in  a  state  as  undisturbed  as 

o 

possible.  Comparing  these  values  with  in  situ  measurements  using  prob 


has  related 


good  agreement  has  been  obtained.  Mackenzie^*” 

®  Hamilton's  measurements  to  observed  bottom  loss.  Despite  the 
great  importance  of  quantitative  knowledge  of  the  ocean  bottom, 
comparatively  little  effort  has  been  devoted  to  this  subject.  Much 
improvement  is  hoped  for  in  the  next  year  or  two.  For  further 
material,  references  30-48  in  the  bibliography  may  be  consulted, 

L  2.  4  Shallow  Water  Propagation 

Mackenzie^’ and  Marsh  and  Schulkin^^  have  made  exten¬ 
sive  analyses  of  the  sound  field  and  its  fluctuations  in  shallow  water. 

9 

Mackenzie  has  developed  a  theory  of  shallow  water  propaga¬ 
tion  based  on  Hamilton  and  Shumway's  description  of  the  bottom  in 
terms  of  density,  sound  speed,  and  absorption.  In  addition,  a  scat¬ 
tering  property  is  superposed  using  a  cos  6  or  co8^0  law  of  scat- 

o 

tering.  Taken  together  with  the  empirical  f  dependence  of  ab- 
0  sorption,  these  expressions  are  used  to  account  for  the  attenuation. 

Sea  surface  scattering  is  neglected.  A  frequency  shift  or  dispersion 
attributed  to  the  bottom  has  also  been  observed  by  Mackenzie  in 
studies  of  fluctuation  statistics  in  shallow  water  propagation. 

Marsh  and  Schulkin  have  described  propagation  in  shallow  water 
through  a  semi-empirical  bottom  loss  factor  which  is  also  dependent 
On  Sea  state.  The  coupling  between  the  bottom  and  the  sea  surface 

arises  because  the  average  angle  of  incidence  of  the  sound  on  the 

o 

bottom  is  a  function  of  sea  state.  This  description  has  the  happy 
faculty  of  being  consistent  with  deep  water  surface  scattering  results 
as  a  function  of  sea  state  and  the  deep  water  bottom  loss  results.  In 
the  first  limiting  ray  zone,  both  the  Mackenzie  and  the  Marsh  and 


8 
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Schulkin  descriptions  are  identical,  despite  the  use  of  the  different 

bottom  loss  curves,  and  the  different  absorption  laws.  There  is 

also  further  theoretical  and  experimental  evidence^^ for  a 

spreading  law  of  r  (cylindrical)  beyond  a  certain  range,  pre- 
-3/ 

ceded  by  a  law  of  r  at  intermediate  ranges,  while  the  spherical 
law  of  r”^  prevails  at  short  ranges.  Williams^^”*^  e.t  al  have  made 
extensive  studies  of  low  frequency  propagation  in  shallow  water, 
also  finding  a  sea  state  dependence.  They  are  now  studying  the 
loss  of  energy  into  the  bottom  through  the  establishment  of  trans¬ 
verse  waves.  Tolstoy^®  has  extended  the  normal  mode  treatment 
of  Pekeris^. 

1.  2.  5  Internal  Fluctuations  in  the  Sound  Field 

The  sound  field  may  be  changed  by  the  phase  interference  of 
existing  rays,  focusing  due  to  the  widening  or  narrowing  of  beam- 
lets,  or  the  addition  or  subtraction  of  new  rays.  Phase  interference 
usually  occurs  when  a  few  rays  following  slightly  different  paths  re¬ 
combine.  If  the  phases  occur  randomly  over  the  interval  of  0  to  2  tt, 
then  the  statistics  of  amplitudes  will  follow  a  Rayleigh  distribution*. 

If  the  signal  is  very  much  larger  than  the  others,  then  a  "Rice"  dis¬ 
tribution^’ (  sine  wave  plus  random  background)  is  obtained.  Changes 
in  focusing  lead  to  normally-distributed  amplitudes. 

The  addition  or  subtraction  of  rays  may  be  brought  about  by 
relative  motion  of  transmitter  or  receiver  beams;  the  motion  of  re¬ 
flecting  surfaces  such  as  the  ocean  surface,  due  to  waves  or  tides, 

e  ® 

or  the  change  in  the  sound  speed  structure  (refraction)  due  to  internal 
waves,  surface  waves,  turbulence,  currents,  or  heating  and  cooling. 
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Such  processes  may  also  cause  phase  changes,  amplitude  changes 
or  even  frequency  shifts  along  single  rays. 

Understanding  fluctuations  in  the  sound  field^°  thus  involves 
an  understanding  of  the  process  of  scattering  (and  reverberation)  in 
relation  to  the  size,  degree  of  discontinuity,  distribution,  and  time 
behavior  of  ocean  inhomogeneities.  In  the  following  paragraphs  we 
summarize  the  situation  with  regard  to  various  types  of  scattering. 

1.  2.  5. 1  Volume  Scattering 

Volume  scattering  contributes  second-order  effects  to  fluctua¬ 
tions  in  field  intensity  and  direction  of  propagation.  Although  a  small 
scattering  attenuation  factor  does  exist,  this  should  become  more  im¬ 
portant  at  great  ranges  and  low  frequencies. 

The  basic  investigations  of  volume  scattering  in  the  ray-theory 
range  were  performed  by  P,  G.  Bergmann^^  who  derived  expressions 
for  the  fluctuations  of  the  acoustic  path  and  the  intensity.  The  wave- 
theory  treatment  of  scattering  was  essentially  formulated  by  C,  L. 
PekeriS*®  and  E.  L.  Carstensen  and  L.  L.  Foldy,^^  who  investigated  the 
effect  of  bubbles  in  water.  D.  Mintzer^^  extended  the  work  of  Pekeris 
and  succeeded  in  deriving  an  expression  for  the  first-order  space- 
average  value  of  the  scattered  intensity.  A.  M.  Obukhov^®  obtained  an 
expression  for  the  fluctuation  of  phase  and  amplitude,  which  in  the  low 
frequency  range  is  equivalent  to  Mintzer's  result  and,  in  the  ray-theory 
range,  is  equivalent  toBergmann's  solution.  Chernov®’  investigated 

O 

the  correlation  distance  for  two  sound  receivers  and  also  wrote  a  book 
on  the  propagation  of  energy  through  random  media.  Other  Russian 

workers  active  in  the  field  are  Krasilnikov,  and  Karauainikou.^®  ’ 

e 


10 


e 


D.  S.  Potter  and  S.  R.  Murphy’^  who  studies  the  effect  of  the  correlation 
functions  e”  ,  and  e  “  in  scattering,  derived  an  expression  for 

the  maximum  phase  delay  of  the  space  average  contributions  of  the 
scattering  patches.  They  also  derived  an  expression  for  the  coefficient 
of  variation  in  tRe  whole  frequency  range. 

Lighthill,  Kraichnan,  Lyon  and  Tatar  ski  have  studied  the 
propagation  of  sound  waves  through  turbulent  media.  Skudrzyk  et 
made  an  important  contribution  to  this  subject  through  their  unified 
theoretical  treatment,  the  introduction  of  the  Kolmogorov  "patch"-size 
distribution  of  temperature  in  the  ocean  and  the  application  of  these 
concepts  to  actual  thermal  and  acoustic  experimental  results.  Skudrzyk 
has  also  contributed  to  the  bearing-error  problem  by  giving  an  expression 
for  the  deviation  in  direction  due  to  focusing  and  interference  fluctuations. 
This  analysis  is  also  applicable  to  volume  reverberation. 

Volume  reverberation,  as  relatej;^  to  the  statistical  distribution 
of  the  scattering  mechanisms,  has  been  treated  by  Carleton,^’  and  also 
Cron  and  Schumacherf®  Lafond'^*’®^  has  been  making  detailed  studies  of 
the  thermal  microstructure.  The  original  work  in  this  area  was  carried 
out  by  Lieberman,®^  and  Urick  and  Searfoss,'^^ 

The  bibliography  of  Skudrzyk  (reference  #77)  may  be  consulted  for 
references  to  the  material  cited  in  this  section. 

1.  2,  5,  2  Internal  Waves 

Lee®®  has  made  computations  of  the  effects  of  internal  waves  on 
sound  field  intensity.  In  addition,  some  study  has  been  made  of  the  propa¬ 
gation  effects  of  the  magnitude  of  the  sound  speed  gradients  in  the  thermo- 

0 

dine.  ® 


1.  2.  5.  3 


Bearing  Error  Studies 


A  satisfactory  theory  of  scattering  should  enable  one  to  study 
the  effects  of  fluctuations  in  horizontal  homogeneity  both  in  bearing 
error  studies,  and  phase  front  fluctuations  and  correlations.  In 
this  connection,  important  work  has  been  done  by  S.  J.  Gershman  and 
G.  L  Tuzhilkin.®*' 

1,  2,  5,  4  Bubbles 

e 

The  existence  and  persistence  of  layers  of  bubbles  arising  from 
the  wakes  of  ships,  and  wind  and  wave  action  have  been  known  to  fur¬ 
nish  an  important  pApagation  effect.  During  World  War  II,  it  was 

found  that  bubble  screens  absorb  and  scatter  non-resonant  sound  waves 
Va 

according  to  an  f  law,  a  dependence  which  has  been  observed  for 
sea  surface  scattering  at  low  grazing  angles.  The  original  treatments 
by  Carstensen  and  Foldy,^^  and  by  Meyer  and  Skudrzyk®^  have  been 
applied  by  Laird  and  Kendig,®®  and  Fox  et  al®’  to  measurements  of  ab¬ 
sorption  and  phase  velocity.  Surface  reverberation  measurements  by 
Urick®^’®®  indicate  that  a  layer  of  bubbles  may  explain  the  observed 
phenomena  of : 

a)  Reverberation  becomes  independent  of  grazing  angle  at 
about  10® 

b)  Pulses  start  returning  to  the  receiver  at  a  depth  of  about 

40  ft.  from  the  surface. 

o  e 

The  subject  of  characteristics  of  bubbles  in  water  continues  to  receive 
increased  attention. 

o  1,  ^  5,  5  Intensity  Computations 

o 

The  use  of  high-speed,  high- storage  computing  machines  has  made 
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ray -tracing  a  practical  and  useful  tool.  There  are  still  some  questions 
raised  by  Pederson?®  concerning  the  propriety  and  necessity  of  rounding 
off  discontinuities  in  sound  velocity  as  applied  to  linear  fitting  of  computed 
sound  speeds  versus  depth.  Intensity  computations  require  the  knowledge 
of  all  rays  which  go  through  two  points  in  addition  to  divergence  compu¬ 
tations  and  diffraction  corrections,  Pederson  has  also  programmed  the 
mode  solution  of  the  wave  equation  as  a  boundary  value  problem. 


Chapter  II 


QUANTITATIVE  DESCRIPTION  OF  THE  SCXJ  ND  FIELD 

2.  1  Introduction 

This  chapter  presents  a  quantitative  method  for  the  calculation 
of  the  average  sound  field.  It  is  based  on  a  simple  representation  of 
the  physical  environment,  tables  and  equations  relating  acoustic  pa- 
rameters  to  the  environment,  and  approxirmte  equations  for  the  sound 
field. 

The  number  of  parameters  which  must  be  stated  in  order  to  des¬ 
cribe  the  average  sound  field  is  rather  large;  this  fact  precludes  the 
feasibility  of  preparing  precomputed  tables  or  nomographs  of  the  sound 
field  itself,  except  possibly  in  special  cases.  It  is,  therefore,  neces¬ 
sary  to  resort  to  a  synthetic  method  for  calculating  the  sound  field.  We 
have  chosen  to  use  ray  geometry  as  the  main  basis  for  this  calculation,  o 
modified  and  supplemented  in  various  special  situations. 

Section  2.  2  contains  a  brief  physical  representation  of  the  envi¬ 
ronment  sufficient  to  permit  calculations  of  the  sound  field.  Equations 
for  the  acoustic  properties  of  this  representation  are  presented  in  Sec¬ 
tion  2.  3  as  functions  of  the  conventional  oceanographic  measures.  In 
Section  2.  3,  detailed  equations  are  given  for  the  sound  field,  and  in  Sec¬ 
tion  2.  4  simple  methods  are  given  for  estimating  the  sound  field  in  the 
special  cases  of  shallow  water,  and  of  sound  channels, 

2,  2  Physical  Environment  of  Underwater  Sound 

2.  2.  1  General  Considerations 

The  physical  laws  describing  ti.e  propagation  of  sound  contain 
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several  functions  which  are  properties  of  the  environment  on  which 
the  sound  is  impressed,  and  others  characteristic  of  the  generation 
and  measurement  of  the  field.  Table  I  presents  these  functions  in 
outline  form.  According  to  this  outline,  the  sound  field,  as  repre¬ 
sented  by  Item  5,  is  determined  as  a  function  of  Items  1  through  4, 
which  in  turn  are  interrelated.  The  outline  states  that  there  are 
more  or  less  well  defined  regions  of  space  which  potentially  influence 
the  sound  field.  In  each  region,  the  complete  time  dependent  equa¬ 
tions  of  state  must  be  known;  these  state  equations  are  divided  into 
physical,  chemical,  and  biological  groups. 

It  is  assumed  that  the  equations  of  state  are  available.  This 
is,  of  course,  not  literally  true.  Qualitatively,  we  know  very  little 
about  suitable  representations  of  the  ocean  floor  and  its  underlying 
materials.  Likewise,  a  similar  statement  can  be  made  about  regions 
possessing  appreciable  ice  covers.  Quantitatively,  we  know  little 
concerning  the  dynamic  character  of  sound  speed  in  the  hydrosphere. 

In  addition,  the  number  of  data  points  necessary  to  specify  a  given 
situation  can  be  very  large,  and  we  are  forced,  for  reasons  of  econo¬ 
my,  to  adopt  simple,  approximate  representations. 

Despite  the  foregoing,  the  results  of  practical  application  have 
showm  that  even  a  gross  representation  of  the  environment  permits  a 
useful  and  reasonable  calculation  of  the  sound  field.  At  this  time,  it 
appears  that  the  important  basic  mechanisms  are  understood  and,  with 
the  exception  of  the  ocean  bottom,  the  necessary  physical  variables 
have  been  fairly  well  determined,  and  the  connections  between  them 
have  been  established. 
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Given  the  necessary  data  corresponding  to  Item  2  of  the  outline, 
the  related  acoustic  properties  {Item  3)  may  be  determined.  These 
are  functions  of  both  the  environment  and  the  generators  of  the  sound 
field  (Item  4.1)  .  The  important  acoustic  parameter  is,  of  course, 
the  frequency.  We  have  admitted  a  dependence  on  wave  amplitude, 
partly  for  generality,  but  also  as  a  reminder  that  the  usual  brand  of 
acoustics  is  an  approximate  linearized  aspect  of  the  larger  subject  of 
hydrodynamics.  Having  made  this  observation,  any  consideration  of 
finite  amplitude  effects  will  be  ignored  in  what  follows.  The  acoustic 
properties;  sound  speed;  absorption;  reflection  and  scattering  coef¬ 
ficients  are  presented  below  in  quantitative  form. 

Having  obtained  the  acoustic  properties,  and  stated  the  sonar 
parameters,  the  sound  field  is  determined  by  the  equations  given  in 
Section  4  of  this  chapter, 

2.  2,  2  Geophysical  Properties  of  the  Environment 

2.  2.  2,  1  Physical  Properties 

These  properties  are  time  dependent  point  functions  in  each 
region.  The  physical  properties  of  velocity,  temperature,  and  pres¬ 
sure  need  no  elaboration.  By  depth,  we  mean  the  function  describing 
the  geometrical  interface  between  the  regions.  As  such,  it  refers  to 
instantaneous  wave  elevation  at  a  geographical  point;  to  the  projection 
of  ice  into  the  water;  to  the  bathymetric  variable  characterizing  the 
ocean  floor;  or  to  substrata  in  the  ocean  bottom.  Of  these,  only  the 
wave  elevation  enters  explicitly  in  this  treatment  of  the  subject. 

For  the  calculation  of  the  mean  sound  field,  the  statistical  proper¬ 
ties  of  the  sea  surface  appear  to  be  sufficient,  because  of  both  theoretical 
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and  experimental  reasons.  A  descriptive  account  of  the  sea  surface 
has  been  given  recently^®  in  connection  with  theories  of  sound  scat¬ 
tering,  The  single  variable  of  significance  is  the  mean  wave  height  h, 
which  is  derivable  from  the  local  wind  speed  for  the  fully  developed 
sea,  or  from  the  wind,  fetch,  and  duration  in  case  of  the  developing 
sea.  In  the  latter  case,  the  methods  of  reference  91,  may  be  em¬ 
ployed  for  determining  both  wave  height  and  the  associated  spectral 
distribution. 

The  influence  of  the  sea  surface  upon  the  average  sound  field 
is  represented  in  terms  of  the  wave  spectrum  A^{to), which  gives  the 
distribution  of  wave  amplitude  among  the  spectral  components  at  angu¬ 
lar  frequency  This  influence  is  described  in  reference  EO.  That 
discussion  is  limited  to  the  fully  developed  sea,  and  we  so  limit  the 
discussion  here.  However,  the  same  methodology  can  be  applied  to 

A  1 

the  co-cumulative  spectra  for  the  non- equilibrium  sea. 

According  to  Neumann  and  Pierson, 

C  2  2^  \ 

A  ^  (w)  =  — r  exp  ■— *  ; — I  cm^  sec  (1) 

w®  -  w‘‘  8^  f 

.  ,  _5 

where  c  =  4,  8  x  lO*  cm  sec  ,  g  is  the  acceleration  of  gravity 
(980  cm/sec^)  ,  s  is  the  wind  speed  in  cm/sec  and  «  the  angular  fre¬ 
quency  in  rad/sec.  In  addition,  the  relation 

h^  =  2,  42  X  10  ^8^  cm^  (2) 

is  employed  for  the  mean  square  wave  height, 

H  =  1.  77h  (3) 

for  the  average  trough  to  crest  wave  height. 
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The  sea  is  presumed  to  be  isotropic,  or  rather,  the  average  scat¬ 
tering  over  all  azimuths  is  considered  in  developing  the  sound  field, 

2,  2,  2,  2  Chemical  Properties 

For  most  practical  purposes,  the  "salinity"  of  sea  water  can  be 
employed  as  the  single  variable  controlling  the  sound  field.  No  chemi¬ 
cal  properties  of  the  atmosphere  or  ocean  bottom  are  known  to  be  sig¬ 
nificant. 

The  salinity,  along  with  temperature  and  pressure,  determines 
both  the  acoustic  phase  velocity  and  absorption  coefficient,  although  the 
various  salts  contributing  to  total  salinity  do  not  contribute  equally  in 
their  affect  upon  these  acoustic  properties.  The  individual  contributions 
are  not  known,  and  as  a  practical  matter,  the  total  salinity  is  employed 
in  determining  the  acoustic  properties.  This  is  reasonably  accurate, 
since  the  ratios  of  the  various  salt  contents  are  remarkably  independent 
of  total  salinityi  except  in  very  brackish  water. 

The  presence  of  air  or  vapor,  either  dissolved  or  suspended  as 
bubbles,  is  known  to  have  an  influence  upon  sound  propagation.  Measure¬ 
ments  of  attenuation  in  water,  "quenched  water"  and  other  related  effects 
were  reported  during  World  War  II,  More  recently,  reflections  from  ap¬ 
parent  bubble  layers  beneath  the  sea  surface  have  been  reported  by  several 
observers.  The  subject  has  not  received  enough  attention,  however,  to 
warrant  including  methods  either  or  predicting  effects  upon  the  sound  field, 
or  for  predicting  the  presence  of  scattering  or  absorbing  regions  of  this 
type  in  the  water.  A  similar  statement  can  be  made  concerning  biological 
effects.  In  considering  other  sonar  properties,  such  as  reverberation, 
the  existence  of  biological  scatter ers  cannot,  of  course,  be  ignored. 


2.  3  Acoustic  Properties 

2.  3.  1  Speed  of  Sound  in  Sea  Water 

2,  3,  1,  1  Introduction 

In  recent  years,  with  the  development  of  an  operational  speed 
of  sound  meter,  a  re-examination  has  been  made  of  the  formulas  de¬ 
veloped  over  the  years  by  less  direct  methods. 

92 

R,  J.  Urick  was  the  first  to  use  an  interferometer  to  make  sound 
speed  measurements  in  the  ocean  during  the  World  War  11,  The  basic 
instrument  which  has  come  to  be  accepted  for  this  purpose  is  a  5  me 
pulse-type  apparatus.  It  measures  the  time  for  a  pulse  to  traverse  a 
fixed  path  between  transmitting  and  receiving  crystals.  This  instrument, 
called  the  "  sing-around"  velocimeter,  was  designed  at  the  National  Bureau 
of  Standards  by  Greenspan  and  Tschiegg.^ 

W,  Wilson^  of  NOL  has  used  this  technique  to  make  the  most  exten¬ 
sive  laboratory  measurements  of  the  speed  of  sound  in  sea  water  as  a 
function  of  temperature,  pressure  and  salinity,  Del  Grosso^^  of  NRL 
made  the  first  set  of  laboratory  measurements  at  atmospheric  pressure 
uncovering  an  error  of  -3  meter  s/sec  in  the  semi- theoretical  tables  of 
Kuwahara,  and  also  in  the  tables  of  Matthews.  Del  Grosso  used  an  inter¬ 
ferometer  technique,  however. 

K.  V,  Mackenzie’'*’^  of  NEL  and  E.  E.  Hays^®  of  WHOI  have  made 
measurements  with  the  NI3S-type  velocimeter  at  sea  as  a  function  of 
depth  and  have  verified  Wilson's  laboratory  measurements  and  resultant 
tables  with  respect  to  depth. 
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These  measurements  point  to  a  new  pressure  effect  term  in 
isothermal  water  of  about  0.0165  m/sec/meter  near  the  surface 
instead  of  0,  0181,  Figures  supplied  by  Mackenzie  show  the  sound 
speed  corrections  as  a  function  of  depth,  between  the  results  of 
Wilson  and  Kuwahara,  They  show  that  these  are  not  linear  and  reach 
a  minimum  at  about  4000  meters  in  <?epth. 

Table  II,  taken  from  Wilson,  summarizes  the  results  at  atmos  - 
pheric  pressure,  for  various  workers  in  the  field. 

Table  II 

COMPARISON  OF  THE  SOUND  FIELDS  OBTAINED 
BY  DIFFERENT  AUTHORS 

T  =  30*C 
P  =  1.033  kg/cm^ 

Author  Date  S  peed  of  Sound(meter8/8ec) 

and  of 

Reference  Determination  Distilled  Water  Sea  Water  S  =  35*/oo 


Matthew 

1939 

1504.  4 

1543, 1 

Kuwahara’® 

1939 

1543.  2 

Del  Grossed® 

1952 

1509. 6 

1546.  2 

Greenspan® 

19  58 

1509.  44 

Wilson^ 

1959 

1509.  66 

1546.  16 

Also  in  the  present  section,  we  reproduce  Wilson's  equation  of 
so’ond  speed  in  sea  water.  We  prefer  this  equation  to  the  equation  sub¬ 
sequently  published  by  Wilson’*,  since  it  represents  the  situation  in 
actual  sea  water  better,  Wilson's  first  equation,  valid  for  salinity  in 
the  range  between  33  and  37  ’/qoi  was  deduced  from  measurements 
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made  on  actual  sea  water  samples.  The  latter  equation  considered 
measurements  made  by  diluting  the  samples  with  distilled  water. 
The  first  equation  has  a  standard  error  of  about  0,  22  m/sec  and  the 
second  a  value  of  0,  30  m/sec. 


Equation  1. 


C  =  1449.  22  tC^.CpiCg.C3^p 


(4) 


where: 

=  4.  6233T  -  5.  4585  x  10  ‘^T^  +  2.  822  x  10"^  T^  -  5.  07  x  10■^T^ 

Cp  0,60518  X  10-^  P  +  1,0  279  X  l0""Pn  3.451  X  10'’ 3.503  X  10"^^P^ 

Cg  =1.391(3-35)  -  7.  8  X  10"MS.35)S  and 

CgTp  =  (3-35)  (-1.  197  X  10-2  T  +  2.  61  X  l0-4p  -  1.96 

X  lO'^'p^  -  2.  09  X  10-*PT)  +  P(.2.  796  x  ’O’^T 

+  1.  3302  X  10“®T^  -  6.  644  x  10‘®T^)  +  P  (-2.  391 

X  lO'^T  +  9.  286  X  10‘‘«t2)  -  i.  745  x  10"1®P^T. 


In  Equation  (4)  ,  P  is  expressed  in  kg/cm^,  T  is  in  ®C,  S  is  in 
parts  per  thousand,  and  C  is  in  m/sec.  The  coefficients  have  been 
rounded  to  give  the  nearest  0,  001  m/sec  in  each  term.  The  rounding 
e  of  coefficients  will  not  affect  the  computation  of  sound  speed  by  mcr  e 
than  0.005  m/sec,  if  the  computation  is  made  within  the  range  of  vari¬ 
ables  considered  here.  To  convert  meters/sec  to  kiloyards /sec,  multi¬ 
ply  by  1,  0936  X  10^ ;  to  convert  meter  s/sec  to  feet/sec  multiply  by  3,2808, 


2,  3, 1,  2  Review  of  Sound  Speed  Computations  from  Specific  Volume  Data 


e  The  differences  between  the  measured  sound  speeds  and  the  pre¬ 
dicted  sound  speeds  at  atmospheric  pressure  were  discussed  by  Del  Grosso’^ 
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and  Beyer'',  To  account  for  the  differences  between  the  computed 
and  measured  sound  speeds  at  other  pressures,  a  review  was  made 
of  the  method  used  for  computing  sound  speeds  from  specific  volume 
data,  Kuwahara  and  Matthews  computed  sound  speeds  in  se^  water 
from  Newton's  formula. 


where  y  is  the  ratio  of  specific  heats,  p  is  the  density,  and  p  is 
the  isothermal  compressibility.  The  density  was  computed  from 
the  formula 


P  =  TT^IP) 

where  p  is  the  mean  compressibility  as  defined  in  Elquation  (8)  below. 
The  true  compressibility,  P  ,  was  found  from  the  mean  compressibility 
by  the  relation, 


^  Y  dP  1  -  Pp  ' 

It  may  be  seen  from  this  equation  that  the  mean  compressibility  is 
defined  by 

n  =  -  i-  (o, 

Vo  (P-Po) 

where  P  =  0,  Substituting  Eqs,  (6)  and  (7)  into  Eq,  (5)  ,  we  obtain 

V  (1-P(x)" 
ftp  P(dp/dP) 
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The  object  in  writing  the  equation  for  sound  speed  in  terms  of  jj. 
instead  of  p  was  to  allow  the  use  of  an  empirical  equation  obtained  by 
V,  Ekman.  In  Ekman' s  equation,  p.  is  written  as  a  function  of  tempera¬ 
ture,  pressure,  and  salinity.  This  empirical  equation  was  based 
principally  on  Amagat’ s  specific  volume  tables  for  distilled  water  and 
a  few  representative  measurements  of  the  specific  volume  of  sea  water 
made  by  Ekman,  The  pressure  dependence  of  p  in  Ekman' s  equation 
relies  upon  Amagat’ s  pressure  measurements.  Therefore,  the  effect 
of  pressure  on  the  speed  of  sound  in  water  should  have  the  same  charac¬ 
teristics  regardless  of  whether  Ekman's  or  Amagat' s  data  is  used.  In 
Equation  (9)  ,  p  can  be  obtained  either  from  Equation  (8)  and  Amagat' s 
data,  or  from  Ekman's  empirical  equation.  If  Equation  (8)  and  Amagat' s 
data  are  used,  the  second  derivative  of  sound  velocity  in  distilled  water 
with  respect  to  pressure  is  positive.  If  Ekman's  equation  for  p  is  used, 
the  second  derivative  is  negative.  It  is  apparent,  therefore,  that  the 
negative  curvature  of  sound  speed  as  a  function  of  pressure  has  been  in¬ 
troduced  by  the  inclusion  of  Ekman's  work.  Since  Ekman's  equation  is 
estimated  to  be  accurate  to  about  three  parts  per  thousand,  sound  speeds 
computed  from  this  equation  cannot  be  more  accurate  than  4,  5  m/sec. 

In  actuality,  the  differences  observed  between  the  measured  and  the  com¬ 
puted  sound  speeds  are  less  than  this  amount.  Comparison  of  the  effect 
of  pressure  on  the  speed  of  sound  is  not  possible  in  other  work  because 
of  the  large  pressure  increments  taken  and  the  small  amount  of  data  ob¬ 
tained  in  the  pressure  range  considered  here, 

2,  3,  1,  3  Note  on  Operational  Experience  With  Sound  V elocimeter 

The  "  sing-around"  velocimeter  is  essentially  an  ultrasonic  delay 
line  which  synchronizes  a  relaxation  oscillator.  The  repetition  frequency 
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is  controlled  by  the  transit  time  in  the  water  for  a  pulse  traveling 
between  two  transducers.  The  transit  time  in  turn  depends  on  the 
speed  of  sound.  High  stability  and  sensitivity  are  feasible  in  the 
system,  permitting  absolute  measurement  of  speed,  to  be  made  with 
an  uncertainty  of  less  than  1  ft/sec.  Readings  taken  in  the  field  are 
reproducible  to  within  about  0,  5  ft/sec. ,  taking  into  account  battery 
aging  and  temperature  effects  on  the  circuit  elements.  Lowering  of 
the  meter  at  the  rate  of  a  fathom  per  second  is  feasible  but  tends  to 
eliminate  detail  in  the  velocity  structure.  Calibration  tests  after  a 
three  month  cruise  have  shown  velocity  measurements  made  with 
the  instrument  to  be  still  correct  within  1  part  in  10,  000.  The  only 
maintenance  required  was  flushing  of  transducers  with  fresh  water 
and  an  occasional  battery  change. 

Comparisons  of  velocimeter  measurements  with  Wilson's  tables 
have  been  reported  by  Mackenzie’^  and  Hay^®.  These  indicate  com¬ 
paratively  large  deviations  near  the  surface  of  the  ocean,  and  offer 
additional  confirmation  to  earlier  evidenc^’ which  has  attributed 
velocity  anomalies  to  the  presence  of  entrapped  air  bubbles.  If  this 
is  the  case,  it  is  still  not  dear  whether  the  observed  effect  would  be 
dispersive,  and  possibly  not  applicable  to  most  sonar  frequencies. 

A  A 

Very  strong  evidence  of  dispersion  has  been  reported  by  Fox  et  al 
in  water  containing  bubble  concentrations  of  2  x  lO'^parts  by  volume. 
The  effect  has  also  been  studied  by  Laird  and  Kendig®*  ,  on  the  basis 
of  calculations  of  phase  velocity  made  from  measured  absorption. 

2.  3.  2  Sound  Absorption  in  Sea  Water 

The  recommended  expression  for  absorption,  a,  in  sea  water 
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as  a  fiinction  of  temperature,  pressure  and  salinity  is  : 

Bt’ 

(l  -  6,  54  X  10  ^  P)  nepers/meter  (10) 


a  = 


SA 


where: 

S  is  the  salinity  in  parts  per  thousand 


f 

A 

B 

P 

T 


is  the  temperature  dependent  relaxation  frequency  in  kc  at 

u  •  Q  in  "152o/(T+  273  )  I 

atmospheric  pressure  =  ZUV  x  10  »•  J 

is  the  acoustic  frequency  in  kc 

is  a  constant  ( 2,  34  x  10”^)  for  the  ionic  relaxation  process 
in  sea  water 

is  nearly  constant  (3,  38  x  10"^)  for  the  pure  water  viscosity 
mechanism, 

is  pressure  in  Kg/cm^  or  atmospheres 
is  the  temperature  in  *C 


To  convert  a  to  db/kyd,  multiply  by  the  factor  7,  943  x  10*, 
2.  3.  2. 1  Discussion 


Apparently  there  is  now  comparatively  little  disagreement  with 
the  temperature  dependence  of  f~.  The  approximate  form  which  it 

^  ^  -D/t 

follows,  and  which  is  plotted  in  Figure  1,  is  f  =  Ce  '  ,  although  the 

D  /r  ^ 

more  accurate  dependence  is  f^  =  ETe  represents  the 

temperature  dependence  of  the  pure  water  viscosity  term.  However, 
it  must  be  emphasized  that  is  not  the  time  constant  associated 

with  the  viscous  process.  Figure  1  compares  measured  values  for  the 
temperature  dependence  of  the  pure  water  relaxation  with  l/f^.  It  may 
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be  seen  that,  except  for  temperatures  above  75°  F,  the  error  in  this 
term  is  less  than  10%,  the  term  itself  usually  being  considerably  less 
than  the  first  terms  for  frequencies  normally  encountered  in  under¬ 
water  applications. 

Murphy  and  Garrison^^  have  made  measurements  over  trans¬ 
mission  paths  in  Puget  Sound  at  frequencies  of  60  kc,  142  kc,  272  kc 
and  467  kc.  They  conclude:  "Although  the  frequency  dependence  of 
^he  absorption  coefficient  given  by  these  measurements  follow  the  nor¬ 
mal  relaxation  law,  the  magnitudes  are  lower  than  those  given  by  Del 
Grosso  by  4  to  9  db/kyd". 

2.  3.  2.  2  Concentration  Dependence 

For  MgS04  solutions  with  concentrations  n  in  range  0,  05  to  0, 10 
mole/liter,  Kurtze  and  Tamm'^  found  a/n  to  be  substantially  constant. 
Their  data,  taken  from  Del  Grosso?^  is  plotted  in  Figure  2,  They  also 
noted  that  sea  water  of  salinity  35  %o  was  equivalent  to  MgS04  solutions 
with  n  =  0.  014.  Leonard  and  Wilsorf”®  investigated  the  concentration 
dependence  for  values  of  n  between  0.  003  and  0,  020,  Their  data  shows 
a/n  to  be  substantially  constant  also,  although  it  is  about  35%  lower. 

Previously,  Beyer^®^  indicated  that  the  excess  absorption  in  MgS04 
solutions  varied  as  the  square  root  of  the  concentration  at  lower  frequen¬ 
cies  and  had  a  "more  complicated  dependence"  at  higher  frequencies.  A 
theory  by  Leontovichf®'*  calls  for  the  absorption  to  be  proportional  to  the 
square  root  of  the  concentration  at  low  frequencies  and  proportional  to 
the  square  of  the  concentration  at  high  frequencies. 
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.  MAXIMUM  ABSORPTION  PER  1IAVELEN67H 


2.  3,  2.  3 


Discrepancy  in  A 


Since  the  linear  dependence  of  absorption  in  concentration  n  is 
in  agreement  for  these  low  concentrations  it  seems  that  the  major 
problem  can  be  traced  to  a  disagreement  in  the  determinaiicn  of  A. 

Thus,  it  appears  that  Leonard  and  Wilson's  laboratory  value 
for  A  is  lower  than  the  laboratory  value  obtained  by  Kurtze  and  Tamm 
for  MgS04  (see  Figure  2) .  Likewise  the  field  values  for  A  found  by 
Murphy  and  Garrison^*  are  definitely  lower  than  those  obtained  by 
Marsh  and  Schulkin^  at  lower  frequencies  in  about  the  same  ratio. 

It  is  felt  that  the  low  saline  waters  of  Puget  Sound  which  Murphy  and 
Garrison  used  could  have  a  lower  percentage  of  MgS04  than  sea  water. 

2.  3.  2.  4  Absorption  Dependence  on  Pressure 

Curiously  enough  the  dependence  of  absorption  on  pressure  in 
sea  water  has  not  been  studied.  However,  studies  have  been  made  of 
the  dependence  of  absorption  in  pure  water  and  also,  associated  and 
non  associated  liquid^°®"^°\  Measurements  on  water  have  been  made 
up  to  pressures  of  2000  kg/  cm^.  Table  III,  reproduced  from  reference 
108,  shows  that  the  absorption  loss  decreases  as  the  pressure  is  raised. 
One  would  expect  similar  behavior  in  sea  water. 

Table  HI 

DEPENDENCE  OF  ABSORPTION  ON  PRESSURE  FOR  WATER,  30"  C 
p(atm)  0  500  1000  1500  2000 

10”  a  f-  18.5  15.4  12.7  1  i.  1  9.9 
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2.  3.  2.  5 


Dispersion 


The  relaxational  absorption  in  sea  water  gives  rise  to  a  change 
of  sound  speed  with  frequency  or  dispersion.  This  effect  is  treated 
in  great  detail  by  Hertzfeld  and  Litovitz?®*  Del  Grosso^^  has  computed 
that  the  maximum  dispersion  in  a  0.  014  molar  solution  of  MgS04 
{equivalent  to  sea  water)  is  about  2-4  crr/sec.  He  used  an  expression 
of  Kneser 


(11) 


Such  an  extremely  small  v^ue  of  dispersion  would  hardly  cause  trouble, 
even  for  systems  with  long  range  communication  applications.  How¬ 
ever,  it  should  be  pointed  out  that  the  dispersion  problem  is  much  more 
severe  than  indicated  above.  For  any  attenuation  curve  that  is  a  function 
of  frequency,  a  dispersion  curve  for  any  four  terminal  network  may  be 
deduced.  Since  the  attenuation  curve  will  include  scattering  losses  at 
the  surface,  losses  at  the  bottom  and  into  the  bottom,  scattering  losses 
in  the  volume  of  the  ocean  leakage  losses  from  sound  channels,  as  well 
as  absorption,  the  problem  can  be  quite  complex  and  important.  Theo¬ 
retically  there  is  a  relation  between  the  temperature,  pressure  and 
salinity  dependence  of  sound  speed  and  absorption  in  sea  water.  If  we 
know  these  relationships  for  one  with  great  precision,  then  the  effect  of 
these  variables  on  the  other  should  be  determined  in  principle. 


2,  3.  3  Acoustic  Properties  of  the  Ocean  Floor 


At  the  time  of  this  writing,  major  attention  is  being  given  to  this 
subject,  and  it  is  not  practicable  to  give  a  useful  accoiant  of  advances  in 
the  general  state  of  knowledge.  It  is,  therefore,  recommended  that  the 
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curves  reproduced  in  Figure  3  be  employed  until  such  a  time  as  they 
can  be  replaced  by  data  giving  dependence  of  the  properties  upon  en¬ 
vironment  and  possibly  other  factors. 

However,  T.  Bell^®’  has  pointed  out  that  the  curves  do  not  take 
explicit  account  of  multipath  effects  in  bottom  reflection,  and  that 
differences  of  as  much  as  six  db  can  exist  depending  upon  the  number 
of  ray  paths  involved  in  a  particular  situation. 


31 


f  9Q)  'SSOI  NOiiOa 


GRAZING  ANGLE.  «  (DEG) 
FIGURES  BOTTOM  VS.  GRAZING  ANGLE 


2.  4  Calculation  of  the  Sound  Field 

2.  4.  1  General  Equations  and  Parameters 

In  this  section,  the  equations  connecting  the  sound  field  with  the 
other  variables  in  the  model  will  be  listed  and  developed  where  neces¬ 
sary. 

The  excess  pressure  P  is  the  variable  to  be  used  in  describing 
the  sound  field,  and  the  fundamental  equation  governing  Pis  the  wave 
equation 

n^P  -  Q'P 

V  at^c.i  (12) 

In  Equation  (12)  ,  ^  ^  is  the  Laplacian  operator  and  c  the  speed 
of  sound.  For  a  unique  determination  of  P,  boundary  conditions  and 
initial  conditions  must  be  stated.  The  boundary  conditions  require  that 
P  and  the  particle  velocity  v  be  everywhere  continuous.  The  initial 
conditions  prescribe  the  distribution  of  some  linear  combination  of  P 
and  V  over  a  closed  surface  (or  set  of  closed  surfaces)  representing 
the  source  of  the  sound  field. 

In  a  fluid,  the  P  and  v  are  connected  by  the  relation  including 
the  density  p: 

Because  of  this,  the  pressure  and  velocity  generally  cannot  be  stated 
independently  for  initial  conditions. 

In  the  following  developments,  the  source  will  be  taken  as  an 
isotropic,  point  source,  with  a  simple  harmonic  angular  frequency  of  oo. 
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In  this  case,  Equation  (12)  becomes 


u> 

where  k  -  — 
c 

the  source, 


and  p  =  Pe 


iwt 


0  (14) 

Near  the  source,  if  R  is  the  distance  to 


P  = 


i(kR  -  wt) 
e 

R 


(15) 


For  other  initial  conditions,  results  may  be  obtained  by  Fourier 
synthesis,  (i.  e.  ,  complex  superposition) . 


The  solution  for  the  sound  field  may  be  obtained  with  the  aid  of 
ray  geometry.  Corresponding  to  Equation  (  12)  ,  there  is  a  set  of  sur¬ 
faces  determined  by  the  equation 

|ye|^  =  (16) 

where  n  -  and  is  the  index  of  refraction  relative  to  the  reference 

c 

speed  Cg. 


The  surfaces,  where  0  is  constant,  constitute  a  set  of  wave  fronts. 
Orthogonal  to  these  wave  fronts  are  skew  curves,  which  are  the  rays  of 
the  sound  field.  If  s  is  the  arc  length  measured  along  a  ray,  the  equa¬ 
tion  of  the  ray  is 


d  ^n  dx  \  0  n 

ds  y  ds  I  9  X 


(17) 


and  similarly  for  y  and  z. 


From  the  manifold  of  rays,  those  rays  connecting  specific  points 
may  be  selected  by  integrating  Equation  ( 17)  and  applying  boundary  con¬ 
ditions.  The  selected  set  of  rays  provide  the  geometry  for  determining 
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the  sound  field  at  the  points  of  interest.  The  equations  for  this  de¬ 
termination  are  presented  in  the  sequeL  There  appears  to  be  no 
point  in  developing  the  subject  in  complete  spatial  generality,  be¬ 
cause  of  the  lack  of  detailed  knowledge  of  the  instantaneous  oceano¬ 
graphic  variables,  and  because  of  the  enormous  complexity  of  any 
calculations  which  might  be  attempted.  Thus,  while  special  situa¬ 
tions  such  as  that  discussed  by  Lee®^  will  continue  to  be  of  interest, 
the  detailed  description  of  the  sound  field  presented  below  will  be 
confined  to  the  case  of  depth  dependent  sound  speed. 

Equation  (12)  is  deriv'able  from  the  hydrodynamical  equations 
together  with  the  equations  of  state  of  the  materials  in  the  environ¬ 
ment,  under  certain  simplifying  assumptions”®.  It  is  basically  valid 
for  static  conditions,  in  which  case  the  parameters  of  the  field  are 
constant  in  time.  This  is  more  particularly  true  of  Equation  ( 14)  . 
However,  if  the  parameters  vary  slowly,  compared  to  the  harmonic 
variation  of  the  source,  then  time  varying  fields  can  be  constructed 
from  a  sequence  of  static  fields.  This  point  will  be  developed  subse¬ 
quently.  Likewise,  spatial  variations  in  the  parameters,  other  than 
in  depth,  can  be  treated  by  connecting  local,  one-dimensional  environ¬ 
ments,  after  the  manner  employed  by  Gardner^” 

The  wave  number  k  is  complex  valued  and  is  generally  not  a 
linear  function  of  the  frequency  as  might  be  implied  by  Equation  ( 14) . 
The  situation  in  sea  water  is  presented  in  Section  2,  3.  1. 

The  boundaries  are  represented  by  the  Rayleigh  reflection  co¬ 
efficients,  and  in  the  case  of  the  ice  free  sea  surface,  by  a  scattering 
coefficient.  Sufficient  information  is  not  available  at  this  time  to  permit 
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the  use  of  more  than  the  reflection  coefficient.  The  bottom  is  to 
be  treated  as  though  it  were  a  fluid,  and  if  velocity  and  attenuation 
data  are  available  for  a  particular  situation,  effects  of  propagation 
into  the  bottom  may  be  included  by  treating  it  as  an  extension  of 
the  hydrosphere. 

2.  4,  2  Ray  Geometry 

There  are  two  methods  of  displaying  the  ray  geometry  pictori- 
all>i  The  first  is  called  the  ray  diagram  and  is  illustrated  in  Figure 
4,  The  ray  diagram  shows  the  paths  taken  by  selected  rayc  Trough 
the  field.  Many  computer  methods  and  programs  have  been  developed 
for  dealing  numerically  with  ray  diagrami^^”,^*^  A  convenient  program 
at  USNUSL,  based  on  the  method  of  Cole,  is  discussed  in  reference 
117. 

For  the  purpose  of  computing  the  sound  field,  a  second  display, 
which  we  called  the  field  diagram,  is  illustrated  in  Figure  5,  Although 
this  class  of  diagram  is  mentioned  by  Horton*^  and  is  widely  used  in 
radio  propagation,  its  application  to  underwater  sound  has  been  limited. 
The  field  diagram  is  employed  to  identify  the  rays  passing  from  the 
source  through  affixed  point  in  the  field.  When  the  rays  have  been  iden¬ 
tified,  the  total  field  may  be  determined  by  combining  the  contributions 
of  each  ray.  Figure  6  shows,  for  the  parameters  indicated  below,  the 
number  of  rays  arriving  at  the  receiver  as  a  function  of  range.  It  is 
clear  that  the  job  of  calculating  the  sound  field  by  this  method  is  not 
trivial.  It  should  be  equally  clear,  however,  that  all  the  rays  contributing 
to  a  point  should  be  considered,  since  the  relative  contributions  of  dif¬ 
ferent  classes  will  vary  from  point  to  point.  Although  many  rays  may 
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FIGURE  4  NON  DIMENSIONAL  RAY  DIAGRAM 
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FIGURE  6.  NUMBER  OF  RAYS  CONNECTING  SOURCE  AND  RECEIVER 
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be  required  per  point,  the  method  is  completely  tractable  in 
comparison  to  other  formulations. 

In  Figure  5,  the  range  between  the  source,  S,  and  the 
receiver,  R,  is  plotted  against  the  ray  inclination  at  the  sur¬ 
face,  for  all  rays  with  the  exception  of  rays  reflected  from  the 
bottom  at  short  ranges.  A  vertical  line  drawn  at  any  range  of 
interest  will  intersect  the  given  curves  at  several  points.  Each 
point  of  intersection  iaentifies  a  particiolar  ray,  providing  the 
ray  inclination  exceeds  the  minimum  value  corresponding  to  the 
horizontal  tangency  at  the  source. 

In  order  to  clarify  Figures  4  through  6,  the  following  ex¬ 
planation  of  the  parameters  used  is  given. 

In  Figure  4,  the  non-dimensional  ray  diagram,  the  hori¬ 
zontal  scale  is  the  range  variable,  u,  which  is  equal  to  gr/2.  The 
vertical  scale  is  the  depth  variable,  v,  equal  to’^Zgz,  In  both 
cases,  g  is  the  fractional  increase  in  sound  velocity,  per  unit  depth. 
The  source  is  located  at  v  =  4/l0  Six  cycles  of  the  ray  tangent 
to  the  source  are  shown  between  the  source  and  the  symmetrical 
receiver  point. 


In  Figure  5,  the  field  diagram,  the  scales  are  equal  to  those 
in  Figure  4.  The  vertical  scale  also  indicates  the  sine  of  the  grazing 
angle  at  the  surface,  (v  =  o)  .  The  heavy  lines  represent  the  source 
depth  (v  =  4/10  and  an  assumed  bottom  (v  =  0.  4)  .  Equal  source 
and  receiver  depths  are  assumed. 
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Figure  6,  number  of  rays  connecting  source  and  receiver, 
is  plotted  with  the  same  parameters  as  Figures  4  and  5,  Bottom 
reflected  rays  are  not  considered  here. 

Temporarily,  take  the  origin  of  cylindrical  coordinates  (r,  z) 
to  be  at  the  source,  and  the  source  inclination  to  be  Gg,  Then  the 

equation  of  a  point  on  the  ray  may  be  written 

z 

cot  fl  d  £,  (18) 

0 

0  is  the  inclination  of  the  ray  at  the  point 

In  Equation  (18)  ,  cot  0  must  be  finite  and  continuous  throughout 
the  interval  of  integration.  In  general,  there  will  be  points  at  which 
one  of  these  conditions  is  not  satisfied.  Furthermore,  the  general 
ray  is  periodic  in  nature,  and  there  A^ill  be  a  set  of  such  points  in 
each  period,  or  cycle  of  the  ray. 

Cot  0  will  be  infinite,  and  the  tangent  to  the  ray  will  be  parallel 
to  the  r  axis  at  points  where  the  sound  speed  is  a  relative  maximum. 

These  points  are  called  turning  points.  The  relative  maxima  can  occur 
at  fixed  depths,  or  at  depths  determined  by  the  parameters  of  the  ray 
in  question. 

Cot  0  will  be  discontinuous  where  a  ray  is  reflected  at  a  boundary. 
Here,  the  depths  of  the  points  of  discontinuity  will  be  fixed.  If  the  bounda¬ 
ries  are  not  perpendicular  to  the  depth  axis,  special  consideration  must 
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be  given  to  the  relationship  between  the  angles  of  incidence,  reflection 
and  transmission.  Several  computer  programs  have  been  developed  to 
account  for  rays  reflected  from  sloping  bottoms  (see  reference  111,  for 
example) .  Williams^^^  gives  an  approximate  method  of  treating  sloping 
bottoms  in  shallow  water.  The  field  diagram  is  particularly  useful  in 
the  ray  solution,  since  one  need  only  draw  in  the  bottom  profile  and 
make  the  proper  adjustments  to  the  rays  affected, 

A  point  of  discontinuity  or  infinity  will  be  called  a  fiducial  point. 
Regardless  of  the  ray  parameters,  there  will  be  exactly  two  depiths 
where  fiducial  points  occur  for  a  given  ray,  and  the  ray  will  cycle  be¬ 
tween  these  depths.  Since  the  origin  of  coordinates  is  taken  at  the 
source  of  rays,  the  origin  will  lie  between  these  two  depths,  including 
the  special  case  where  the  origin  is  at  one  of  the  depths.  These  depths 
will  be  denoted  Z+  and  Z,  ,  It  is  obvious,  then,  that  Z+  *  Z  $  Z., 

Let 

F^(z)  =  I  coted^  (19) 

Here  e  may  be  +  or 

r  =  2F+  (0)  +  2F.  (0)  (20) 

has  been  termed  the  "cycle  range"  by  Cole,“^  and  is  twice  the  skip- 
distance  defined  by  Marsh  and  SchulkinV 

Equation  (18)  can  now  be  written 

r=r,+  [F+(0)  -  F^.  (z)]  +  rj  _^.(0)  -  F_  (z)  ]  +  2mF^.  (0)  +  2nF.  (0) 

(21) 

The  integers  m  and  n  are  the  numbers  of  contacts  with  the  fiducial 
points  Z-t-  and  Z_  respectively,  and  the  integers  rjg  take  values  1,  0,  -1 
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according  to  the  following  table. 


0  5^  0 

0  0 

>1+ 

1 

-1 

Z  ^  0 

0 

0 

Z  «  0 

0 

M 

n_ 

-1 

■■ 

The  smallest  of  the  integers  m,  n  is  the  number  of  complete  cycles 
in  the  ray.  If  Gq  >o,  m  ^  n  and  vice  versa.  In  any  case,  jm  -  nj  <  1. 


In  addition  to  the  purely  geometrical  properties  of  the  rays  the 
eikonal  K  is  important.  We  have 


K 

u 


/  T  +  iA 

0  ^ 


(22) 


T  is  the  "Travel  time,"  and  A  represents  absorption.  Corresponding 
to  the  quantity  Fg  ,  define 


=  I 


dt. 


sin  9  cos  0 


’  H  = 


f 


sin  0  cos  0 


;23) 


Then, 


CqT 

cos  00 


Cq  a 
cos  6  0 


=  r}+  Jg+(0)  -  G+(z)j  +  n.  |^G_(0)  -  G_(z) 
+  2mG  ^  (0)  +  Zr]  G_  (0) 


r7+  ^  (0)  -  H+  (z)  j  +  I^H-' 

+  2m  H^(0)  +  2r}  H_(0) 


(0)  -  H_(z) 
(24) 
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It  is  of  interest  to  note  that,  by  virtue  of  the  mean  value  theorem 


of  the  calculus, 


T 


r 

c • cos  6 


(25) 


where  c  and  cos  0  correspond  to  some  point  along  the  ray.  This  equation 
shows  the  extreme  possible  departure  of  the  travel  time  from  that  which 
would  be  associated  with  the  corresponding  refraction  free  situation. 

2,  4.  3  Primary  Field 

Any  point  in  the  sound  field  reached  by  one  or  more  rays  will  be 
called  a  primary  point,  and  the  field  produced  at  these  points  by  the  rays 
will  be  called  the  primary  field.  The  total  field  is  the  sum  of  the  primary 
field  and  the  field  due  to  scattering  and  diffraction.  The  primary  field 
will  almost  always  be  dominant  at  primary  points,  as  discussed  in  Section 
2.  4.  6.  At  other  points,  the  secondary  effects  must  be  considered,  and 
are  treated  below  in  Sections  2.  4,  4  and  2,  4.  5, 

Take  the  origin  of  coordinates  at  the  mean  level  of  the  sea  surface, 
with  the  source  on  the  polar  axis,  and  let  the  cylindrical  coordinates  of 
the  source  be  (0,0,  Zg) ;  that  of  a  field  point  be(r,  4>,  Z)  as  in  Figure  7. 

Let  Fj  be  any  ray  connecting  the  two  points  and  0i  the  inclination 
of  the  ray  at  the  source.  The  pressure  associated  with  Fj  will  be  desig¬ 
nated  Pi. 

The  total  primary  field  is  the  complex  sum 

p  =  2i  Pi 

over  all  rays  reaching  the  field  point.  The  equation  for  Pi  will  be  stated. 


44 


and  then  the  individual  terms  discussed. 

Pi  =  qJ(5l  )  S“  (26) 

Each  term  in  this  equation  is  a  function  of  the  ray  parameters 
and  so  a  subscript  (i)  is  implicit. 

The  quantity  q  represents  the  optical  approximation  to  the  pres¬ 
sure,  and  is  equivalent  to  the  pressure  reduction  due  to  geometrical 
divergence  of  the  ray  bundle  as  discussed  by  Horton^*,  among  many 
others. 
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The  quantities  S  and  B  are  the  Rayleigh  reflection  coefficients 
for  the  ocean  surface  and  floor,  respectively,  m  the  number  of  surface 
contacts  and  r)  the  n\amber  of  bottom  contacts.  If  the  ray  penetrates 
either  surface  and  returns  by  refraction  or  further  reflection,  the  re¬ 
flection  coefficients  must  be  replaced  by  the  proper  combination  of  re¬ 
flection  and  transmission  coefficients  in  an  obvious  manner. 

The  quantity  S  is  sufficiently  close  to  -1  for  nearly  all  purposes, 
and  in  the  sequel  this  will  be  assvimed.  The  state  of  knowledge  con¬ 
cerning  B  has  been  presented  in  Section  2.  3.  3.  B  is  a  function  of  the 
grazing  angle  of  the  ray  at  the  bottom,  the  acoustic  frequency,  and 
possibly  the  environment.  The  fiinctional  dependence  is  largely  em¬ 
pirical  at  this  time,  although  Mackenzie^®  has  shown  that  the  acoustic 
properties  of  sediments  may  be  reduced  to  first  principles  when  suf¬ 
ficiently  detailed  knowledge  of  the  material*  is  available.  We  write 

B  =  B  (u,  Gg.c) ,  (31) 

8b  being  the  grazing  angle  at  the  bottom  and  e  indicating  potential  en¬ 
vironmental  dependency, 

n.  is  the  specular  scattering  coefficient  of  the  sea  surface^®.  For  small 

scattering,  3A  iAq 

n  =  1  -  0.  485  b  H  sin  Gg  (32) 

where  H  is  wave  height  in  feet,  b  is  equal  to  w  H/27r,  and  8g  is  the 

grazing  angle  at  the  sea  surface.  We  take  the  same  parametric  depen  - 

dence  to  hold  for  large  scattering,  so  that,  in  general, 

a  =  n  (b^^^  H  ^^‘’sinGg)  (33) 

Figure  8  shows  measured  values  from  many  sources, with  Equation  (32) 
plotted  as  a  solid  curve  in  the  form,  surface  loss  =  -  10  log  n. 
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ffiEOUEICY  X  SAVE  HEIGHT  F.H.  FT.  X.C. 
FIGURE  8.  SEA  SURFACE  SCATTERIRG  LOSS 


The  dotted  curve  is  the  recommended  extension  of  Equation  (33) 
to  use. 

This  curve  has  been  drawn,  taking  into  account  proper  weighting, 
and  the  additional  fact  that  the  scattering  should  level  off  at  high  fre¬ 
quencies. 

The  irreversible  conversion  of  acoustic  to  thermal  energy,  or 

-A 

absorption,  is  represented  by  the  term  e  ,  in  which 


2.  4.  4  Scattered  Field 

Sound  scattering,  as  a  propagation  effect,  produces  three  basic 
effects.  These  include; 

a)  contribution  to  the  primary  field 

b)  creation  of  sound  energy  in  shadows 

c)  contribution  to  fluctuation 

The  latter  effect  is  discussed  in  Chapter  3,  The  other  effects 
are  primarily  caused  by  the  ocean  surface  and  bottom.  Although  the 
bottom  has  scattering  properties,  they  are  not  sufficiently  understood 
at  this  time  to  permit  quantitative  calculations.  The  scattering  proper¬ 
ties  of  the  sea  surface  have  been  discussed  by  Skudrzyk'^^,  among  others, 
and  the  model  employed  by  Marsh^^®  et  al  will  be  employed  here  for  esti¬ 
mating  the  scattered  field. 

The  scattered  field  is  analyzed  into  its  plane  wave  components, 
and  for  an  incident  plane  wave.  In  terms  of  the  geometry  of  Figure  9, 
Marsh  gives,  for  small  scattering. 
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r 

where : 


d-v^') 

I 


2/2  -  - 
0.  098  H  ®  (  1  -  V*  )  exp  ( -  1  /aM  ^ 

6.  33  H  b  (  1  -  vM  ^ 
2-(v-y)^-2  cos  ijj 
r  d  Xd  p 


(35) 


as  the  ratio  of  scattered  intensity  in  the  direction  (X,  p)  to  incident 
intensity.  For  application  to  a  definite  calculation,  integration  over 
(X,  p)  must  be  effected.  This  is  prohibitively  difficult.  It  would 
appear,  however,  that  the  case  of  greatest  interest  is  that  in  which 
propagation  is  supported  by  repeated  surface  reflection.  In  that 
case,  the  calculation  can  be  materially  simplified  by  some  approxi¬ 
mations.  We  neglect  multiple  scattering  and  joint  refracted- scat¬ 
tered  processes.  Consider  Figure  10, 

We  calcvilate  the  scattered  field  at  the  mid-point  of  the  cycle 
and  use  it  as  an  average  over  the  entire  cycle.  Since  multiple  scat¬ 
tering  is  neglected,  only  primary  rays  need  be  considered. 


The  scattering  process  entails  no  loss  of  energy,  but  rather  a 
redistribution  in  direction  of  propagation.  Thus,  the  total  scattered 
intensity  will  equal  that  lost  from  the  specularly  reflected  ray.  Ac¬ 
cordingly,  the  scattered  field  due  to  the  primary  ray  in  question  is 
approximately 

,  CFe’^^ 

—  (36) 


R‘= 


where 


C 


sin  p  d  p 


e 
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The  results^  give  empirical  expressions  for  the  sound  field 
below  the  isothermal  layer,  which  should  be  equivalent  on  the  average 
to  those  obtained  from  Equation  (36)  ,  where  they  can  be  compared. 

The  data  is  limited  to  frequencies  above  two  kilocycles  per  second. 
Below  this,  the  scattered  field  will  be  rather  small.  Furthermore, 
at  about  this  frequency,  the  effects  of  diffraction  start  to  become  im¬ 
portant. 

2.  4,  5  Diffraction  Field 

The  diffraction  features  of  the  sound  field  at  points  reached  by 
rays  has  been  considered  above.  At  other  points  (shadows)  the  dif¬ 
fraction  field  can  be  predominant. 

Shadows  always  have  a  boundary  which  is  a  limiting  ray  having 
a  turning  point  as  a  sound  speed  maximum.  In  the  shadow,  the  dif¬ 
fraction  field  will  be  nearly  constant  along  lines  parallel  to  the  limiting 
ray,  and  will  diminish  exponentially  in  proceeding  away  from  the  limi¬ 
ting  ray.  In  Figure  11  there  are  two  limiting  rays,  and  there  will  be  a 
diffraction  field  associated  with  each. 


The  development  leading  to  Equation  (29)  is  applicable  in  prin¬ 
ciple  to  the  present  situation.  However,  no  useful  results  can  be  ob¬ 
tained,  In  case  the  sound  speed  is  a  linear  function  of  depth,  and  if 


y  ^ 


_1 

c  dZ 


(37) 


then  the  diffraction  field  pressure  is  nearly  equal  to 

1/3 

-3(tt) 

Pi  « 

H  and  V  are  the  horizontal  and  vertical  distances  shown  in  Figure 
p^  is  the  pressure  at  the  limiting  ray. 


(38) 


1 1  and 
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There  appear  to  be  few  instances  that  Equation  (38)  can  be 
applied  with  reliable  results.  Reference  117  indicates  that  the 
shadow  zone  field  is  much  larger  than  that  calculated  from  the 
above  equation.  On  the  other  hand,  the  equation  is  surely  appli¬ 
cable  if  the  prevailing  assumptions  are  justified.  Two  definite 
possibilities  exist: 

1)  There  are  other  sources  of  sound  in  the  shadow  zone, 

2)  The  sound  speed  structure  is  not  accurately  known. 

These  can  be  discussed  for  the  cases  shown  in  Figure  11, 

When  the  turning  point  is  at  the  surface,  and  if  bottom  reflections 
do  not  contribute  appreciably,  the  diffraction  field  should  dominate. 
Morse*^*  has  shown  the  qualitative  effect  to  be  expected  from  small 
changes  in  the  speed  gradient  near  the  surface.  This,  coupled  with 
recent  evidence  in  near  surface  sound  speed"  indicates  that  a 
reliable  calculation  of  the  field  can  be  anticipated  if  the  speed  struc¬ 
ture  is  definitely  known. 

When  the  turning  point  is  at  depth,  as  under  a  surface  sound 
channel,  the  scattered  field  can  be  expected  to  predominate  at  the 
higher  frequencies  where  measurements  exist.  At  the  lower  frequen¬ 
cies,  a  reasonable  calculation  of  the  diffracted  field  can  be  expected 
using  Equation  (  38 )  , 

2.  4.  6  Normal  modes 

At  sufficiently  low  frequencies,  the  methods  of  calculation  based 
on  ray  geometry  are  no  longer  valid.  The  use  of  normal  modes,  or 
solutions  based  upon  separation  of  the  wave  equation  can  be  utilized  in 
such  cases.  As  a  practical  matter,  the  normal  mode  solution  will  be 
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particularly  useful,  and  will  differ  appreciably  from  the  ray  solution, 
when  there  is  a  single  mode  which  is  propagated  without  substantial 
diffraction  losses. 

The  basic  theory  and  many  references  to  the  normal  mode  theory 
are  presented  by  Officer  in  reference  122,  Since  the  appearance  of 
that  publication,  some  work  has  been  published^^’  indicating 

generally  good  agreement  between  theoretical  and  experimental  results. 

Other  than  the  shallow  water  mode  field,  the  only  application  of 
practical  interest  appears  to  be  to  surface  sound  channels.  However, 
at  low  frequencies,  surface  sound  channels  are  seldom  important  in 
determining  the  total  field.  On  this  account,  there  does  not  seem  to  be 
sufficient  interest  to  develop  general  equations  and  procedures  for  the 
application  of  normal  mode  soliitions. 

As  an  aid,  however,  the  results  of  Voorhis***  are  reproduced  in 
Figure  12,  This  shows  the  net  attenuation  due  to  diffraction  of  the  first 
mode  for  various  values  of  the  parameters.  When  this  attenuation  is 
appreciable  compared  with  other  attenuations  (scattering,  reflection  and 
absorption) ,  then  consideration  should  be  given  to  the  use  of  normal 
mode  calculations.  For  that  purpose,  Pedersen^^*  has  made  certain 
corrections  and  improvements  to  the  work  of  Marsh,  His  resvilts  are 
quoted  below,  since  they  represent  the  only  known  complete  results 
necessary  for  the  calculation  of  the  sound  field  itself,  rather  than  the 
features  of  it,  such  as  are  presented  in  the  other  works  cited. 

The  computation  is  based  on  four  physical  parameters  and  three 
parameters  determined  by  the  experimental  set-up.  The  physical 
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FIGURE  12.  ATTENUATION  DUE  TO 
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10  000 


FIRST  MODE 


parameters  are:  Cg ,  the  sound  velocity  at  the  sea  surface;  Zj,  the 
depth  of  the  surface  channel:  and  yg  and  \i,  the  gradients  in  and  below 
the  channel  respectively.  In  this  study  vg  is  positive  and  yi  is  negative. 

The  other  three  parameters  are  the  frequency,  f,  and  the  source  and 
receiver  depths.  In  this  theory  the  source  and  receiver  depths  are  in¬ 
terchangeable  —  i,  e.  reciprocity  holds.  The  shallower  of  these  two 
depths  is  designated  by  Z^,  and  the  deeper  by  Z^. 

Certain  useful  combinations  of  these  parameters  are: 

ta  =  tb  -  ^b/^ii  -  Vo/Vi 

1/3 

r=(jrfyo")  /cg,  M  =  2cg  ZjF  Vyo 
where  N  =  greatest  integer  <  (2M  ri)  +  .  25  ,  the  number  of  trapped 

modes. 

When  Marsh's  equations**®  are  evaluated  in  detail  the  following  form 
results  for  the  product  of  the  depth  functions: 

(ta)-  U  (tfe)  =(M^j  ).  F(ta,tt  )/D  (39) 

D  will  be  discussed  later,  F(ta,t^)  takes  one  of  three  possible  forms  as 
follows: 

Case  I  -  when  both  source  and  receiver  are  in  the  channel;  i.  e,  t^  ^  t|j  <  1 

A{ta)*  A(tb)  (40) 

Case  n  -  when  both  source  and  receiver  are  below  the  channel;  i.e.  l^ta^tb 

Case  ni  -  when  one  is  in  and  one  below  the  channel;  i.  e.  t^  ^  1  ^  tb 

A  (ta)  h2(B(;^)  (42) 
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where 


A(t)  =  hi  (MX^  )  hz  (MX^  -  Mt)  -  hi  (MX^)  hi  {MX^  -  Mt)  (43) 

The  hi  and  h^  fvinctions  are  solutions  of  Stokes'  equation  and  are  discussed 
at  length  and  tabulated  in  reference  119.  First  derivatives  of  these  func¬ 
tions  are  designated  as  hi^  and  h^',  is  the  complex  eigen  value  and 

will  be  given  presently. 

Bt  -  pMMXn  -  M)  +  IM  (44) 

P 

A  (1)  and  Bi  are  obtained  by  setting  t  =  1  in  the  above  expressions. 

D=(p^-1)  [b^  +{MX,,  -  M)  AMd]  -  2.12429  30  (  45) 

where  the  negative  constant  is  the  square  of  the  Wronskian, 

W  (hih^)  and  B  is  defined  by  B  =  hj  {MXr))h2'(MX,,  -  M)-  h2{MX,,)hitMX^  M) 

(46) 

MXq  is  a  complex  root  of  the  characteristic  equation 

ph2(Bi).  B-  h/  (Bi).  A(l)  =  0  (47) 

In  the  case  of  strong  trapping  the  following  approximation  to  the  roots 
may  be  useful. 

Re  (MXr,)=[-^i— (n  -  \ 

Im(MX,^)=  I  |^Re(MX^)j"  exp  ^^l  .  pr’)  (M  - ReMX^)* 

(48) 

The  velocity  potential  ((>  is  given  by 

iwT  ^ 

<!>  =  -  TT  ie  r)  U^(ta)-U^(tb)  (49) 

Af,  =  (Tn  -  i  T  Tt 

0^  ~  K  -  r  Re  (MXi,) 
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c 


Trj  ~  rim(MX^) 


Co 


The  damping  coefficient  is  given  by  T  p  while  0”^  is  needed  in  the 
phase  determination. 

For  sufficiently  large  values  of  r 

Ho^(X^  r)  ~  (--jf-T-)  exp  i  (  (Tn  r  -  -  f^r"^  (50) 

The  transmission  loss  H  is  given  by 
H  =  -  10  log  1  <i>  p 

From  Equations  (49  )  and  (50) 

fr  N  T  2  Tzl 

1<1>)  ^  ^  j  2  cos  0^  +  S  sin  0^^  V/r  (51) 

where 

=  |u^(ta)  U^(tb)l  ®  (52) 

0T7  =  -  (Tn-r 

i|;q  =  arg*  (tb)} 

2.  5  Practical  Calculation  of  the  Sound  Field 

2.  5.  1  General  Considerations 

For  detailed  and  precise  calculation,  the  methods  of  Section  2.  4, 
or  others  of  equivalent  scope,  appear  to  be  ue  cessary.  There  are  two 
practical  difficulties,  however,  which  demand  that  simpler  methods  be 
employed  when  circumstances  warrant.  In  the  first  place,  since  the  en¬ 
vironmental  model  is  idealized,  there  is  often  a  lack  of  sufficient  environ¬ 
mental  data  to  permit  complete  employment  of  the  model.  Secondly,  one 
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is  often  more  interested  in  the  average  sound  field,  or  in  general 
aspects  of  the  field  which  can  be  stated  approximately.  This  is 
especially  so  in  systems  engineering,  and  in  strategic  problems, 
where  only  the  gross  structure  of  the  sound  field  can  possibly  be 
considered. 

Two  situations  exist  which  have  great  practical  importance, 
and  which  can  be  characterized  rather  simply  with  fair  approxima¬ 
tion.  These  are  sound  channels  and  shallow  water.  Both  are  de¬ 
fined  below,  and  methods  given  for  the  approximate  determination 
of  the  sound  field. 

2.  5.  2  Sound  Channels 

A  sound  channel  is  defined  to  be  the  region  of  space  in  which 
the  cyclic  character  of  the  rays  is  predominant,  and  in  which  the 
deepest  fiducial  point  is  a  turning  point.  In  other  words,  the  range 
is  several  multiples  of  the  cycle  ranges,  and  the  influence  of  the 
ocean  bottom  is  ignored. 

Equation  (21)  shows  that,  when  the  source  and  receiver  depths 
are  equal,  the  range  is 

r  =  2mF^  (o)  +  2nF  _  (o) 

(53 

=  mro  +  2  (n  -  m  )  F.(o) 


Hence  when  the  range  is  an  integral  multiple  of  the  cycle  range, 

r  =  mro 


(54) 


dr  _  dr^) 
d  00  "  ""  d  8o 
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If  a  ^  o  we  have  a  relative  maximum  in  the  sound  speed;  if  a  =  oo 
we  have  a  cusp  in  the  speed  profile.  The  case  o  <a<^  can  include 
relative  maxima  if  the  slope  of  the  speed-depth  curve  is  discontinuous 
at  the  turning  point.  The  case  a  =  •«  is  of  no  particular  interest  and 
will  be  disregarded. 


We  are  concerned  with  the  behavior  of  the  derivatives  of  F  near 
a  turning  point.  Using  higher  order  terms  as  z  — >»•  Z  ,  we  may  derive 
the  relationship 


sin  0  — — 

de 


Bin 


(59) 


as 


o,  which  yields 


q 


i( 

e 


ojt  -  kr) 


2  r  ‘sin  0of 


(60) 


From  Equation  ( 60)  we  obtain 

qJ(9)- 


j/  i  ( w  t  -  kr  +  *^12), 

/2  _  /  1.  i 

0 


2jr)  "  e  /  k  n  a 


,V6 


/3(3) 


r^  sin  00- 


cos  80 


We  note  in  passing  that  if  =  1  -  y  Z,  the  case  discussed  by 
Marsh'^^  when  2  a  =  ^  , 


and 


ro 


2  sin  00  cos  Op 
cos^  Op 


(61) 

(62) 


which  is  a  well  known  result 
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The  preceding  discussion  applies  to  isolated  points  in  the  Sound 
channel.  The  general  situation  can  be  estimated  as  follows.  Firstly, 
we  observe  that  both  of  the  quantities  q  and  J(9  )  tend  toward  values 
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o 


representing  spherical  spreading  along  any  given  ray,  as  the  distance 
along  the  ray  increases.  This  is  trivially  true  for  non-cyclic  rays, 
or  for  reflected,  non-refracted  rays.  For  cyclic  rays  the  result  is 
true  because,  as  r  increases. 


r  =  mrn  +  e 


d  r  ^ 

15;=  “-Se;'' 


jLl 

de„' 


m 


d^ro 

dBo" 


(63) 


+  e? 


in  which  e,  Cj,  and  have  maximum  values,  independent  of  r. 

Hence  for  points  chosen  at  random  along  a  ray,  the  probability  that 
j  r  -  mro  |  <  E  tends  to  1  as  r  increases,  where  E  is  any  preassigned 
number.  Equivalent  statements  hold  for  dr /d0o  and  d^r/ilGo^.  We  can 
thus  say  qualitatively  that  every  ray  is  ultimately  spherical.  Now,  in  a 
sound  channel,  the  average  number  of  rays  reaching  a  field  point  in¬ 
creases  as  the  range  increases.  Hence  the  total  field  ultimately  follows 
a  cylindrical  spreading  law.  More  precisely,  if  N  is  the  number  of  rays, 
then 


N 


(64) 


The  quantity  N  can  be  determined  by  constructing  a  field  diagram 
(Figure  5)  ,  or  can  be  estimated  as  follows.  Let  tq  and  ro^  ^e  the  maxi¬ 
mum  and  minimum  cycle  ranges  for  all  possible  angles.  Then 


2.  5.  3 


1  t 


«  1  + 


Shallow- Water  Transmission^^ 


(65) 


Accounts  of  shallow-water  transmission  have  recently  appeared  in 
the  literature.  In  particular,  references  2,  55,  58,  122,  128  discuss  the 
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boundary-valueproblems  associated  with  the  normal  mode  solutions, 
while  references  49,  50,  51,  129  discuss  solutions  at  higher  frequen¬ 
cies  involving  ray  approximations. 

The  principal  problem  in  calculating  the  sound  field  in  shallow 
water  is  the  lack  of  detailed  knowledge  of  the  shallow  water  environ¬ 
ment.  However,  highly  variable  sound  speed  is  the  rule,  and  this 
has  an  important  bearing  on  high-frequency  transmission.  The  proper¬ 
ties  of  the  bottom,  important  at  any  frequency,  are  controlling  at  the 
lower  frequencies. 

Even  given  the  detailed  parameters  of  a  particular  environment, 
the  variability  to  be  expected  is  such  that  calculation  will  be  very  com 
plex  and  tedious.  There  is  thus  the  need  for  comparatively  simple 
equations  representing  the  average  sound  field,  while  retaining  depen¬ 
dence  upon  the  principal  features  of  the  environment. 

A  normal  mode  treatment  is  necessary  if  the  ratio  of  water  depth 
to  wavelength  is  not  large.  In  this  case,  the  accounts  of  Williams**'  *** 
cited  above  appear  to  represent  the  state  of  knowledge  rather  well,  and 
hence  this  specific  aspect  of  the  problem  will  not  be  treated  further  here. 

If  the  water  depth  exceeds  about  four  wavelengths,  then  Mackenzie's, 
along  with  other  results  (not  published)  may  be  expressed  in  a  relatively 
compact  semiempirical  way  compatible  with  a  general  theory  of  under¬ 
water  sound  propagation.  Reasonable  assumptions  are  made  about  the 
mean  character  of  the  shallow-water  environment. 

Both  the  results  of  Mackenzie  and  others  indicate  no  systematic,  de¬ 
pendence  upon  the  depth  of  source  or  receiver.  Propagation  at  extended 
range  is  supported  by  repeated  bottom  and  surface  reflection,  regardless 
of  the  thermal  conditions.  Thus,  there  is  a  strong  surface-bottom  coupling 
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such  that  the  propagation  losses  are  controlled  by  the  number  of  con¬ 
tacts  of  rays  with  both  surface.  The  surface  scatters  the  rays,  where¬ 
as  the  bottom  absorbs  them.  The  thermal  structure  of  the  water  affects 
propagation  through  its  influence  on  skip  distance  and  the  number  of  sur¬ 
face  and  bottom  contacts.  The  propagation  loss  is  thus  represented  in 
terms  of  sea  state  (wave  height)  bottom  type  (or  bottom  loss,  if  known) 

water  depth,  and  isothermal-layer  depth. 

1 

A  skip  distance,  H  =  +  DJ  ^  kyd  is  defined  for  water  depth 

D  an<f  layer  depth  L  in  feet.  If  the  range  R  between  source  and  receiver 
is  less  than  H,  then  the  propagation  loss  is 

N  =  20  log  R  +  QRi-60-kL  db;  (66) 

k^  1  the  near  field  anomaly,  represents  the  mean  contribution  to  the 
field  of  the  multiple  bottom  and  surface  reflections.  It  is  a  function  of 
bottom  loss  and  surface  loss  such  that 

kL  =  10  log  ( 1  +  2r  ^r^  +  r^  +  tk  )  db,  (67) 

and  being  the  surface  and  bottom  intensity- reflection  coeffici¬ 
ents,  respectively.  The  absorption  coefficient  in  sea  water  (in  kyd)  is 
a  and  is  discussed  in  references  9  and  14. 

N  =  1 5  log  R  +  aR  +  a-j-  (  R/H  -  1 )  +  5  log  H  +  60  -  kj^  db.  ( 68 ) 

For  long  ranges,  R  3^  8  H, 

N  =  lO  logR  +  aR  +  a,j,(R/H  -  l)  +  lO  logH  f  64,  5  -  k^db. 

(69) 

These  equations  represent  the  gradual  transition  from  spherical 
spreading  in  the  near  field  to  cylindrical  spreading  in  the  far  field. 

The  loss  coefficient  (db  per  cycle  of  bottom  and  surface  reflec¬ 
tion)  is  represented  as  follows; 
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Table  IV,  Shallow-water  attenuation  (db /bounce) 


Sea 

state  0  1  2  3  4  5 

f 

kcps  Sand  Mud  Sand  Mud  Sand  Mud  Sand  Mud  Sand  Mud  Sand  Mud 


0. 

1 

1. 

0 

1.  3 

1.  0 

1. 

3 

1.0 

1.  3 

1. 

0 

1. 

3 

1. 

0 

1. 

3 

1. 

0 

1. 

3 

0. 

2 

1. 

3 

1.  7 

1.  3 

1. 

7 

1.  3 

1.  7 

1. 

3 

1. 

7 

1. 

3 

1. 

7 

1. 

4 

1. 

7 

0. 

4 

1. 

6 

2.  2 

1.  6 

2. 

2 

1.  6 

2.  2 

6 

2. 

2 

1. 

7 

2. 

4 

2. 

2 

3. 

0 

0. 

8 

1. 

8 

2.  5 

1.  8 

2. 

5 

1.9 

2.  6 

2. 

2 

3. 

0 

2. 

4 

3. 

8 

2. 

9 

4_ 

0 

L 

0 

1. 

8 

2.  7 

1.9 

2. 

7 

2.  1 

2.9 

2. 

6 

3. 

7 

2. 

9 

4. 

1 

3. 

1 

4. 

3 

2. 

0 

2. 

0 

3.  0 

2.  4 

3. 

5 

3.  1 

4.  4 

3. 

3 

4. 

7 

3. 

5 

5. 

0 

3. 

7 

5. 

2 

4.  0 

2. 

3 

3,  6 

3.  5 

5. 

2 

3.  7 

5.  5 

3. 

9 

5. 

8 

4. 

1 

6. 

2 

4. 

3 

6. 

4 

9. 

0 

3. 

6 

5.  3 

4.  3 

6. 

3 

4.  5 

6.  7 

4. 

7 

6. 

9 

5. 

0 

7. 

3 

5. 

1 

7. 

5 

10. 

0 

4. 

0 

5.9 

4.  5 

6. 

8 

4.  8 

7.  2 

5. 

0 

7. 

5 

5. 

2 

7. 

8 

5. 

3 

8. 

0 

Table  V.  Near -field  anomaly  k  (db) 


Sea 

State  0  1  2  3  4  5 

f 


kcps 

Sand 

Mud 

Sand 

Mud 

Sand 

Mud 

Sand 

Mud 

Sand 

Mud  Sand 

Mud 

0.  1 

7.  0 

6.  2 

7.0 

6.  2 

7.0 

6.  2 

7.0 

6.  2 

7.0 

6.  2  7.  0 

6.  2 

0.  2 

6.  2 

6.  1 

6.  2 

6.  1 

6.  2 

6.  1 

6.  2 

6.  1 
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In  the  case  where  upward  refraction  occurs,  it  is  expected  that 
the  sea  state  would  have  a  strong  influence  on  the  propagation.  Actu¬ 
ally,  the  loss  per  contact  looks  like  the  bottom  loss,  rather  than  like 
the  sea- surface  loss  in  that  the  surface- scattering  loss  has  a  much 
larger  frequency  dependence.  This  behavior  with  sea  state  might  be 
expected  if  most  of  the  energy  scattered  from  the  surface  is  reflected 
by  the  bottom  and  returned  for  propagation  down  the  channel.  In  deep¬ 
water,  Og  measures  the  energy  which  is  scattered  out  of  the  surface 
channel  and  which  is  either  lost  in  absorption  on  the  way  to  the  bottom 
or  absorbed  in  the  bottom  through  high  grazing-angle  incidence. 

The  problem  then  is  to  determine  the  extent  of  the  coupling  be¬ 
tween  the  surface  and  the  bottom.  If  rg  is  the  surface-reflection  co¬ 
efficient,  then  ( 1  -  rs)  is  the  surface- scattering  coefficient  and  Ug  = 

-  lO  logiotg  is  the  surface  loss  in  db/bounce.  The  surface  compo¬ 
nents  Tg  and  (1  -  Tg)  must  suffer  different  interactions  with  the  bot¬ 
tom  if  any  sea- state  dependence  is  exhibited.  We  would  expect  scat¬ 
tered  rays  to  suffer  greater  loss  because  they  are  steeper.  The  sim¬ 
plest  expression  satisfying  this  requirement  and  the  observed  data  was 

^t  =  ^a  1  )  ’^b'  (70) 

where  rj.  is  the  fraction  of  energy  transmitted  down  the  channel  when 
a  bottom  event  is  coupled  with  each  surface  reflection,  =  -  lOlogjQr). 
{  shallow- water  attenuation  constant,  db/bounce)  ,  r^  is  the  bottom  re¬ 
flection  coefficient  for  the  angles  of  incidence  occurring  in  shallow 
water  and  is  given  in  Table  IV  for  sand  and  for  mud  bottoms,  -  10 

^  b ^  * 

The  observed  behavior  for  Uj  has  profound  significance  in  the  in¬ 
terpretation  of  the  shallow-water  propagation  mechanism.  The  fact 
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that  the  grazing-angle  rays  represented  by  rg  must  be  multiplied  by 
r^j  indicates  that  there  is  a  bottom  loss  suffered  by  near  grazing  rays. 

The  fact  that  the  scattered  rays  must  be  multiplied  by  r^^  means  that 
the  angular  dependence  of  bottom  reflection  loss  is  such  that  these 
steeper  rays  suffer  twice  the  loss  of  the  near  grazing  rays. 

The  observed  bottom  loss  at  grazing  angles  seems  to  imply  a 
mode  -changing  process  that  takes  energy  for  a  retracted-ray  system 
and  converts  it  to  a  simple  bottom  and  surface  bounce  mode  at  great 
ranges.  This  is  equivalent  to  the  normal  mode  treatment  of  propaga¬ 
tion  in  shallow  water.  This  mode-changing  process  seems  to  occur 
for  any  sound- speed  strvicture.  The  sound- speed  structure  determines 
the  skip-distance  zones.  It  would  appear  that  after  the  first  few  zones, 
the  skip-distance  has  a  purely  local  significance  representing  a  limiting 
free  path  for  rays  which  happened  to  be  scattered  into  grazing  angles 
at  that  range.  If,  as  is  shown  in  the  data,  grazing  rays  must  suffer  one 
bottom  contact  loss  per  skip  zone,  then  they  cannot  be  continuous  from 
one  zone  to  another.  The  possibility  that  horizontal  rays  suffer  a  bottom 
loss  is  hardly  conceivable.  The  transition  from  one  propagation  mode 
to  the  other  probably  occurs  over  several  skip  zones  and  we  have  adopted 
the  convenient  model  that  it  occurs  over  the  first  eight  skip  zones.  When 
the  foregoing  expressions  are  used  to  compute  propagation  loss  for  all 
the  shallow  water  data,  representing  about  100,  000  observations,  then 
Table  VI  gives  the  probable  errors  as  a  function  of  range  and  frequency. 
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Chapter  III 


OTHER  EFFECTS 

o 


3.  1  Introduction 

The  requirements  of  modern  and  probable  future  sonar  systems 
are  such  that  the  variability  of  the  oceanographic  parameters  in  time 
and  space  can  become  an  integral  part  of  the  sound  field  problem. 
Fluctuations  of  the  sound  field  in  the  ocean  are  intimately  related  to 
the  changing  structure  of  the  transmission  path  or  paths  between  the 
transmitter  and  receiver.  How  we  treat  these  fluctuations  in  intensity, 
phase,  direction  resolution,  frequency,  or  pulse  width  depends  entirely 
on  what  result  is  to  be  accomplished. 

If,  for  example,  we  are  trying  to  communicate  over  long  range 
paths,  in  the  design  of  an  optimum  system,  we  must  consider  (a)  the 
time  and  space  variability  of  signal  levels  for  the  modes  of  propagation 
which  can  exist,  (b)  multiple  arrivals  and  multipath  propagation  which 
arises  from  the  physical  nature  of  the  problem,  and  (c)  frequency  sen¬ 
sitive  effects  such  as  dispersion,  absorption,  and  time  and  space  cor¬ 
relation  of  signal  levels,  * 

3,  2  Mean  Variability  of  Transmission  Loss 

The  practical  application  of  the  formulas  for  transmission  loss 
in  Chapter  2  usually  requires  an  estimate  of  the  associated  variability. 
This  is  necessary  in  specifying  performance  measures  of  equipment, 
and  in  optimizing  system  design,  especially  with  regard  to  the  selection 
of  the  best  operating  frequency. 

While  work  is  in  progress,  notably  that  of  Skudrzyl^^,  Mintzer^^, 
and  Tatar  ski^^’in  the  ocean  medium,  with  the  fluctuations  in  the  acoustic 
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parameters,  we  are  forced  to  accept  the  enviroirment  as  it  occurs. 

Thus,  specification  of  the  variability  in  transmission  loss  can  be 
expressed  using  empirical  data  taken  from  experiments  made  with 
the  type  of  system  to  be  studied. 

For  example,  in  previous  attempts  to  correlate  acoustic  be¬ 
havior  with  the  oceanographic  variables,  the  variables  were  con¬ 
sidered  to  be  quasistatic  over  the  acoustic  paths.  The  propagation 
behavior  was  described  in  terms  of  mean  values.  The  propagation 
modes  considered  were  those  which  minimized  the  number  of  classes 
necessary  to  account  for  the  modes  occurring  most  of  the  time  in 
most  locations.  However,  there  remained  a  residual  variability  in 
propagation  loss  within  each  class.  No  particular  study  was  made 
relating  the  variability  within  a  class  with  the  variability  of  the  oceano¬ 
graphic  parameters.  This  is  a  perfectly  feasible  approach  for  the  pre¬ 
diction  of  design  information  for  the  conditions  of  the  experiments. 

If  the  sound  field  is  computed  by  the  methods  of  Chapter  2,  the 
mean  variability  can  be  determined  from  probable  errors  of  the  acous¬ 
tic  properties.  Figure  13>  shows  the  probable  error  in  the  bottom  loss. 
TPk  probable  error  of  the  surface  loss  may  be  determined  from  Figure  8. 

In  addition  to  boundary  variability,  there  is  variation  due  to  the 
body  of  water  itself.  This  is  evident  even  at  short  range,  and  can  be 
estimated  using  the  methods  of  Skudrzyk,  As  a  practical  matter,  the 
inherent  volume  variability  is  about  3  db  in  signal  level. 

If  the  joint  statistics  of  the  i.iultiple  reflection  processes  were 
known,  the  variance  of  the  composite  sound  field  could  be  calculated. 
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FREQUENCY  (KC) 

FICURE  13  PROBABLE  ERROR  VS.  FREQUENCY 


In  the  absence  of  this  knowledge,  two  extreme  calculations  can  be 
made.  One  assumes  independence  from,  and  the  other  assumes 
complete  dependence  upon  the  probable  errors  associated  with  each 
reflection  process.  If  the  variances  computed  in  this  manner  are 
0” ^  and  re  spectively,  then  the  actual  variance  should 

be  bounded  IJ  ^  (J*  ^  <  0^^. 

An  expression  for  the  variance  of  the  propagation  loss  in 
db  is:  1 

(j”=  lO  log  (ratio  equivalent  of  boundary  los8)^j  db  (7l) 

This  overestimates  the  situation,  but  is  a  proper  estimate  for  conser¬ 
vative  calculations. 

At  the  other  extreme,  we  have  simply  ^ 

2 

(f  =  [3^  +  (Tg  +  (Tg^]  (72) 

as  an  example  of  the  minimum  variance  of  the  propagation  loss  for  a 
field  dominated  by  surface  and  bottom  reflected  paths,  and  ^  g^i 
being  the  variances  associated  with  single  refractions  at  surface  and 
bottom,  respectively.  Some  tables  giving  ftie  gross  variability  of  the 
sound  field  have  been  published.  Table  VI  gives  the  probable  error 
of  propagation  loss  for  the  average  shallow  water  environment  along  the 
east  coast  of  the  United  States.  Mackenzi^'*  has  reported  results  for 
shallow  water,  and  has  commented  on  the  nature  of  the  distribution  function. 

3.  3  Factors  Influencing  Variability 

If  the  parameters  of  the  sound  field  were  constant  in  time,  there 
would  be  no  variability.  Changes  in  the  geometry  of  the  field,  such  as 
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motion  of  source  or  receiver,  or  changes  in  the  environment  are 
responsible  for  the  changes  in  the  sound  field  from  its  mean  value. 

Calculations  of  the  fluctuation  to  be  expected  from  thermal 
variations  range  from  microscopic  effects,  through  intermediate 
effects  of  internal  waves  to  the  seasonal  influences  which  are 
readily  observed.  There  is  a  continuum  of  periods  of  environmental 
changes  ranging  from  seconds,  or  fractions  of  seconds  through  sea¬ 
sons  or  even  years.  These  changes  are  entirely  in  the  body  of  the 
hydrosphere  and  its  surface. 

The  effect  of  the  environmental  changes  is  imposed  upon  the 
sound  field  through  changes  in  the  elementary  conotituents  of  the 
field.  In  general,  the  composite  field  is  a  superposition  of  contribu¬ 
tions  from  a  number  of  energy  paths,  which  pass  through  different 
volumes  of  water  and  interact  separately  with  the  sea  surface.  It  is 
the  changing  interaction  of  these  paths  which  causes  variability. 

An  obvious  example  is  that  of  a  ship  suspended  hydrophone, 
moving  about  as  the  ship  moves  in  the  seaway.  Such  motion  can 
cause  enormous  variation  in  the  sound  field.  Even  with  source  and 
receiver  fixed  in  sjiace,  there  is  still  variabilityl^*^  Measurements 
show  a  mean  variation  of  some  2db  in  the  rms  pressure,  and  indicate 
a  more  or  less  well  defined  period  on  ilie  order  of  7  seconds.  This 
variability  is  most  likely  attributable  to  surface  waves. 

In  shallow  water,  with  source  and  receiver  fixed  above  the  bot¬ 
tom,  tidal  periods  have  been  observed. 

The  static  sound  field  ir  the  presence  of  an  instantaneous  internal 
wave  pattern  has  been  discussed  by  Lee^f  A  potential  basis  now  exists 
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for  estimating  the  variability  of  the  sound  field  to  be  expected  be¬ 
cause  of  internal  waves. 

Both  the  characterization  and  the  measurement  of  variability 
are  very  difficult  matters  to  deal  with.  Still,  it  is  only  because  of 
variability  that  the  modern  interest  in  space-time  coherence  of  the 
sound  field  has  developed.  It  appears  possible  to  calculate  varia¬ 
tions  in  the  sound  field  from  a  sequence  of  static  field  calculations, 
in  which  the  environment,  as  represented  by  the  model,  is  varied 
in  accordance  with  the  physical  processes  giving  rise  to  the  environ¬ 
mental  changes.  These  processes  are  of  two  basic  types:  a)  flux 
of  energy  into  the  environment,  and  b)  relaxation  of  energy  concen¬ 
trations  within  the  environment  (increase  of  entropy) ,  In  the  first 
case,  energy  flows  into  the  sea  through  the  surface,  because  of  wind 
action  and  solar  heating,  and  through  the  bottom  because  of  thermal 
gradients  in  the  structure  of  the  earth.  In  addition,  energy  is  supplied 
by  body  forces  (gravitational),  represented  in  the  tides  and  in  Coriolis 
effects  represented  in  current  patterns.  Secondly,  the  relaxation  ef¬ 
fects  represent  a  redistribution  of  energy  stored  from  earlier  eras, 
and  are,  on  a  time  scale,  very  long  compared  with  the  periods  of  varia¬ 
tion  due  to  energy  flux. 

It  would  appear  that  the  predominant  effects  giving  rise  to  varia¬ 
bility  are  wind,  the  sea  surface  heating  (cooling)  and  tides.  Steady 
currents  do  not  produce  variability.  The  time  scale  of  environmental 
changes  associated  with  these  effects  range  from  seconds  to  hours  to 
days  to  seasons.  When  a  periodic  change  in  the  environment  is  impressed 
upon  the  sound  field,  there  generally  will  be  many  periods  produced  in 
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the  sound  field  variation,  and  these  will  be  sub-multiples  of  the  driving 
period. 

Two  predominant  effects  are  present.  One  may  be  the  changing 
phase  relationships  between  different  paths  contributing  to  the  sound 
field.  If  so,  the  fundamental  period  of  fluctuation  can  be  a  very  small 
fraction  of  the  driving  period  when  many  complete  cycles  of  phase  are 
varied  during  the  environmental  cycle.  If  the  phase  change  is  one  cycle 
or  less,  the  fundamental  period  will  equal  the  driving  period.  The  other 
type  of  effect  occurs  when  changes  in  the  shape  or  slope  of  the  reflecting 
or  refracting  surface  introduces  a  varying  convergence  or  divergence. 

In  this  case,  the  driving  period  usually  will  equal  the  period  of  the  acous¬ 
tic  field. 

3.  4  Other  Measures  of  the  Sound  Field 

There  are  many  properties  of  the  sound  field  which  are  not  con¬ 
veniently  or  practically  related  to  the  mean  field  and  its  variability. 

Some  of  these  are  discussed  below.  The  subject  will  be  more  generally 
covered  in  Chapter  IV,  under  new  research. 

3.  4,  1  Impulse  Response  Function 

In  some  communications  applicatiow ,  we  are  interested  in  matching 
the  transmitted  and  received  signal.  Dispersion  properties  of  the  medium 
could  seriously  affect  the  performance  of  such  a  correlation  system.  One 
way  of  obtaining  this  information  is  to  interpret  the  impulse  response  of 
the  propagation  path  as  a  four -terminal  network,  or  filter.  This  will  fur¬ 
nish  phase  and  attenuation  information  as  a  function  of  frequency  for  the 
path.  The  change  of  the  impulse  response  with  time  and  oceanographic 
parameters,  would  then  be  studied  and  charted  for  the  oceans  much  like 
the  acoustic  properties  of  the  bottom.  Pet  erson^^has  shown  that  Horton's 
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attenuation  formula  can  be  ufled  in  this  fashion,  to  account  for  the 
impulse  response  of  the  bottom  to  bomb  shots. 

3.  4,  2  Space  and  Time  Correlation  Functions 

For  some  transmission  problems,  the  size  of  the  receiving 
array  is  limited  by  the  space-correlation  function  of  the  wave  front, 
i,  e.  how  far  aj>art  can  we  put  hydrophone  elements  and  still  have 
the  signals  received  in  phase  ?  Time  correlation  of  signals  at  the 
same  receiver  is  also  important  in  signal  processing, 

3,  4.  3  Concurrent  vs  Sequential  Statistics 

In  some  systems,  a  single  pulse  may  be  transmitted  involving 
several  discrete  frequencies.  In  this  case,  we  would  be  interested 
in  the  correlation  of  signals  from  one  frequency  to  another  at  the 
same  instant  (concurrent  statistics)  and  from  one  instant  to  another 
at  the  same  frequency  (  sequential  statistics), 

3.  4,  4  Multipath  Interference 

The  tendency  of  the  ocean  to  confuse  signals  by  combination  of 
the  images  being  reflected  from  the  surface  and  bottom  will  require 
study  also.  Thorp  and  Power have  shown  the  nature  of  the  problem 
including  the  effect  of  mixed  modes,  or  leakage  arrivals.  There  are 
two  effects  to  be  considered  here: 

3.  4.  4.  1  Effective  Noise  Level  Due  to  Multipaths  Arrivals  (Coherent) 

If  four  equal  reflected  and/or  refracted  rays  arrive  simultaneously 
with  the  direct  signal,  as  might  happen  for  low  frequencies  over  long 
paths,  then  the  signal  interference  could  be  as  much  as  6db  above  the 
signal. 
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3.  4,  2  Effective  Pulse  Lengthening  (Incoherent) 

The  effective  pulse  lengthening  can  be  taken  as  the  time  difference 
between  the  first  arrival  and  the  last  arrival  which  is  less  than  20  db  be¬ 
low  the  signal. 
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Chapter  IV 


REQUIREMENTS  FOR  FUTURE  RESEARCH 
4.  1  General  Comments 

One  way  to  appraise  requirements  for  future  research  is 
to  examine  the  material  presented  in  the  earlier  sections  of  this 
document.  We  shall  do  this,  employing  the  outline  offered  in 
Table  I,  as  a  method  of  organizing  the  subject.  In  so  doing,  we 
shall  try  to  present  priorities,  based  both  on  need  in  terms  of  the 
ultimate  consumers,  and  on  the  scientific  adequacy  of  present  in¬ 
formation. 

One  may  ask  first,  "Is  Table  I  an  acceptable  outline  ?  "  We 
believe  that  it  is,  at  the  present  time,  and  that  future  developments 
will  at  most  add  items  to  the  outline,  without  disturbing  its  struc¬ 
ture.  This  does  not  imply  that  certain  sections  might  be  completely 
reorganized.  For  example,  the  geophysical  properties  might  in  the 
long  run  be  described  statistically  rather  than  deterministically. 

We  will,  in  any  event,  regard  the  outline  as  acceptable,  and  base 

research  requirements  on  the  need  for  providing  information  to  fill 

up  the  data  points  implicit  in  the  outline.  o 

4.  2  Regions 

This  item  of  the  outline  is  simply  a  list.  As  written,  volumes  and 
boundaries  are  listed  as  semantically  equivalent.  There  is  some  logi¬ 
cal  objection  to  this;  a  more  serious  objective  is  that  the  sea  surface, 
for  example,  is  not  a  mathematical  boundary,  but  a  region,  more  or 
less  well  defined,  which  may  include  subsidiary  regions  such  as  the 
boundary  layer  topside,  and  turbulent  air-water  regions  bottomside. 
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Equivalent  statements  can  be  made  for  the  ocean  floor,  and  for  ice 


cover. 


Another  aspect  of  the  regional  division  concerns  turbulence 
paths  within  the  hydrosphere,  as  well  as  "  scattering  layer s".  Are 
these  to  be  regarded  as  identifiable  regions,  with  peculiar  proper¬ 
ties,  or  are  they  automatically  represented  if  the  state  variable  of 
the  ocean  are  given  ? 

Questions  of  this  type  are  not  entirely  a  matter  of  taste,  but 
do  in  fact  depend  upon  the  physical  measures  employed  in  describing 
the  regions.  Furthermore,  commonness  of  nomenclature  amonq  the 
various  disciplines  which  concern  themselves  with  the  different 
regions  is  important  in  developing  the  unified  work  necessary  for 
understanding  the  cross  discipline  of  underwater  sound. 

We  feel,  therefore,  that  thought  devoted  to  the  organization  of 
physical  data  on  the  earth  environment  can  be  fruitful,  and  that  it  is 
an  interdisciplinary  problem. 

4.  3  Geophysical  Properties 

This  section  of  the  outline  refers  to  the  "meat"  of  the  data. 
Herein  falls  the  quantitative  description  of  the  environment.  The 
following  appear  to  be  the  most  important  measurements  needed: 

Depth  Detailed  and  descriptive  profiles  of  the 

ocean  bottom,  and  of  ice  cover. 

®  Velocity  Current  data  at  depth. 

Temperature  Precision  measurements  in  thermal  mi- 

crostructures. 

Chemical  Properties  The  chemical  (and  related  physical) 

properties  of  the  water  just  below  the 
sea  surface,  A  redetermination  of  the 
salt  contents  of  various  sea  waters. 

Biological  Properties  Habits,  distribution  and  physical  proper¬ 
ties  of  all  forms  of  marine  biology. 
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4,  4  Acoustic  Properties 

Given  the  geophysical  properties,  we  coirapute  the  acoustic 
properties.  The  overwhelmingly  important  lactc  today  is  con¬ 
cerned  with  the  effects  of  the  ocean  bottom  upor*  sound  propagation. 

We  can  say  that  the  bottom  reflection  and  scattering  coefficient 
should  be  determined,  but  what  we  mean  is  that  vays  must  be  found 
to  describe  the  bottom  sufficiently  to  predict  is  effects. 

The  second  requirement  is  to  determine  ttie  dependence  of 
sovind  speed  upon  factors  other  than  temperature,  pressure  and 
salinity.  This  should  include  near  surface  effects  as  well  as  posj 
sible  associated  dispersion.  In  connection  with  such  a  study,  it 
would  be  highly  desirable  to  relate  the  coefficients  (temperature,etc, ) 
of  sound  speed  and  absorption. 

4.  5  Parameters 

Again,  this  is  simply  a  list.  Possibly,  ho-wever,  there  are 
other  parameters  more  suitable  than  these,  Peirhaps  the  impulse 
function  is  more  suitable  than  the  frequency  function.  We  will  con¬ 
sider  this  type  of  question  below. 

4,  6  Sonar  Properties 

The  great  emphasis  today  is  on  detail.  Th-e  words  heard  most 
often  are  "space-time  correlation", "multiple  orier  probability  distri¬ 
bution",  and  "fluctuation",  among  others.  Very  often,  one  gathers 
that  the  need  is  for  much  more  "data".  It  is  important  to  recognize 
that  there  is  a  definite  trend  in  progress,  from  tine  circumstances 
where  average  propagation  loss  was  a  decisive  fa-ctor  in  acoustic  prob¬ 
lem,  to  the  situation  where  much  more  must  be  fcnown  about  the  sound 
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field.  It  is  even  more  important  to  recognize  that  this  trend 
is  qualitative,  in  the  sense  that  one  wishes  to  know  new  kinds 
of  things,  more  than  it  is  quantitative,  in  the  sense  that  more 
quantitative  data  of  familiar  form  is  required. 

If  more  detailed  understanding  of  the  sound  field  is  re¬ 
quired,  the  measures  employed  are  going  to  have  to  conform  to 
some  kind  of  conservation  laws,  which  can  provide  powerful  re¬ 
lationships  vastly  simplifying  the  description  and  interpretation 
of  the  field. 

Conservation  of  energy  is  perhaps  the  most  important  single 
principle  in  physics,  and  it  is  curious  that  it  has  found  so  little 
application  in  underwater  acoustics.  Thus,  although  energy  units 
are  frequently  employed  in  the  reduction  of  transient  phenomena, 
this  is  but  a  formal  convenience.  On  the  other  hand,  Uricl^**  has 
shown  how  conservation  of  energy  leads  to  a  natural  extension  of 
the  sonar  equations  to  cover  transient  signals.  One  difficulty  is, 
of  course,  that  the  acoustic  energy  or  intensity  do  not  obey  simple 
differential  equations.  However,  there  may  be  related  quantities, 
such  as  discussed  by  Wolf^^^which  are  conserved  and  obey  a  wave¬ 
like  equation. 

In  addition  to  conservation, it  would  be  highly  desirable  for 
measures  to  be  subdivided  as  little  as  possible.  A  considerable 
loss  of  valuable  interpretation  can  arise  because  of  artificial  division 
of  the  field  into  "reverberation",  "  scattering", 'transmission",  "tar¬ 
get  strength",  and  the  like,  even  though  these  terms  have  a  certain 
qualitative  value.  The  single  function  "Impulse  Response"  has  a 
unifying  significance  in  describing  the  total  result  of  applying  a  driving 
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force  to  the  ocean. 


4.  7  Summary 

In  summary,  the  following  research,  in  order  of  importance, 
is  recommended, 

(1)  Find  new  and  pragmatic  measures  of  the  sound  field. 

(2)  Find  suitable  methods  for  describing  the  acoustic  ef¬ 
fects  of  the  ocean  bottom  and  establish  these  effects, 

(3)  Obtain  detailed  profiles  of  the  ocean  bottom  and  of  ice 
cover. 

(4)  Investigate  the  chemical  (and  related  physical)  proper¬ 
ties  of  water  just  below  the  sea  surface, 

(5)  Determine  the  dependence  of  sound  speed  upon  factors 
other  than  temperature,  pressure  and  salinity, 

(6)  Make  a  redeterrnination  of  the  salt  contents  of  various 
sea  waters, 

(7)  Establish  the  habits,  distribution  and  physical  proper¬ 
ties  of  marine  biological  forms. 

(8)  Study  the  organization  and  methodology  associated  with 
geophysical  measurement. 
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